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ABSTRACT 

The  U.S.  Army  Wartime  Ammunition  Distribution  System  (WADS)  will  experience 
an  unprecedented  demand  for  ammunition  under  the  operational  concept  of  Airland 
Battle.  To  meet  demand,  proper  storage  facility  location  and  an  efficient  flow  through 
the  distribution  network  will  be  required. 

Using  information  from  Army  Field  Manuals,  maps  and  simulation  data  for  de- 
mand, both  a  mixed  integer  program  (MIP)  and  a  sequential,  optimization-based 
heuristic  are  developed  to  model  the  WADS.  The  Generalized  Algebraic  Modelling 
System  is  used  to  implement  both  models.  The  sequential  heuristic  locates  ammunition 
facilities  with  a  binary  integer  program  and  then  directs  ammunition  through  those  fa- 
cilities utilizing  a  network  flow  model  with  side  constraints.  The  MIP  integrates  location 
and  flow  decisions  in  the  same  model.  For  solving  a  typical  scenario,  involving  the  lo- 
cation of  21  storage  facilities  and  the  allocation  of  flows  for  30  time  periods,  the  se- 
quential heuristic  took  22  CPU  seconds  on  an  IBM  3033AP  mainframe  computer.  For 
locating  the  same  number  of  facilities  but  allocating  flows  for  only  3  time  periods,  the 
MIP  took  87  CPU  seconds.  The  heuristic  solution  was  always  within  2%  of  optimality 
on  all  test  problems  that  were  small  enough  to  solve  with  the  MIP. 

Results  indicate  shortcomings  in  the  WADS  as  it  currently  exists.  The  models  and 
analysis  show  that  current  doctrine  is  infeasible  unless  there  is  an  increase  in  lift  assets 
at  the  Corps  level  storage  facilities  and  a  reduction  in  inventory  goals  at  the  Ammunition 
Supply  Points. 
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THESIS  DISCLAIMER 


The  reader  is  cautioned  that  computer  programs  developed  in  this  research  may  not 
have  been  exercised  for  all  cases  of  interest.  While  every  effort  has  been  made,  within  the 
time  available,  to  ensure  that  the  programs  are  free  of  computational  and  logic  errors, 
they  cannot  be  considered  validated.  Any  application  of  these  programs  without  addi- 
tional verification  is  at  the  risk  of  the  user. 
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I.      INTRODUCTION 

The  U.S.  Army  Wartime  Ammunition  Distribution  System  (WADS)  is  presently 
based  on  heuristics  or  "rules  of  thumb"  developed  from  the  experience  gained  during 
World  War  II,  Korea  and  Vietnam.  The  concept  of  Airland  Battle,  developed  since 
Vietnam,  has  generated  concern  within  the  Army's  ammunition  analytical  community 
because,  using  this  concept,  a  greater  demand  for  ammunition  will  be  placed  on  the 
WADS. 

This  thesis  develops  an  integrated  model  which  concurrently  locates  ammunition 
facilities  and  determines  proper  ammunition  flow  using  a  mixed  integer  program.  How- 
ever, because  of  limited  solution  capabilities,  a  sequential,  optimization-based  heuristic 
is  developed  which  decomposes  the  problem.  The  sequential  heuristic  is  composed  of 
two  optimizing  submodels,  an  ammunition  storage  facility  location  submodel  and  an 
ammunition  network  flow  submodel.  Both  models  are  developed  to  test  current  doctrine 
and  provide  a  tool  for  analysis  of  future  systems. 

A.     AIRLAND  BATTLE 

The  concept  of  Airland  Battle  is  simply  a  means  to  defeat  a  large  armored  force 
which  attacks  by  echelon,  through  a  narrow  breach  sector  in  an  opponent's  front;  Soviet 
bloc  forces  are  the  attackers  and  the  Western  forces  are  the  defenders. 

The  idea  of  Airland  Battle  is  to  defeat  the  armored  force's  first  two  echelons  at  the 
proposed  breach  point  and  then  use  a  deep  attack  to  disrupt  any  further  flow  of  the  third 
and  successive  echelons  forward.  Deep  attack  might  be  USAF  fighter  aircraft  or 
bombers,  conventional  cruise  missiles,  special  operations  forces,  tactical  nuclear  weap- 
ons, etc.  This  concept  differs  from  prior  ideas  since  our  forward  units  are  now  required 
to  defeat  two  echelons  rather  than  the  one  echelon  of  past  conflicts.  Defeating  two 
echelons  will  require  more  ammunition  per  unit  engaged  with  the  enemy.  Ammunition 
consumption  will  increase  (Figure  1  on  page  2). 

The  concept  of  Airland  Battle  is  dynamic  since  it  requires  our  widely  spread  forces 
to  move  to  a  breach  point  and  concentrate  assets  to  defeat  the  first  two  enemy  echelons. 
Firepower  far  superior  to  that  used  in  past  conflicts  will  be  required  to  defeat  the  enemy 
at  the  breach  point.  To  achieve  superior  firepower,  a  greater  demand  for  ammunition 
will  be  required  than  in  past  wars.  Ammunition  must  arrive  forward  in  the  type  and 
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Figure  1.      Airland  Battle 


quantity  required;  a  stockout  situation  during  an  enemy  attack  would  most  likely  be 
fatal. 

The  ammunition  community  claims  that  even  in  past  conflicts, 

"Ammunition  availability  constrains  combat  power  before  shortages  of  combat  ve- 
hicles, crews,  maintenance,  repair  parts,  and  POL  [Ref.  1J." 

So,  concisely,  the  problem  is  to  move  ammunition  forward  to  meet  demand  in  the  most 
expeditious  manner  possible  subject  to  available  assets  and  subject  to  constraints  on 
vulnerability  which  result  from  large  inventory  build-ups. 

B.  PROBLEM  SCOPE  AND  PURPOSE 

This  thesis  examines  the  problems  associated  with  locating  ammunition  storage  fa- 
cilities and  determining  an  efficient  flow  of  ammunition  between  storage  facilities  within 
a  generic  Corps  in  a  Theater  of  Operations  (TO).  Network  structure  above  the  Corps 
level,  i.e.  Port,  Theater  Storage  Area,  and  associated  transportation  links,  will  not  be 
considered  due  to  program  and  computer  size  restrictions.  In  addition,  Port  and  TSA 
facilities  are  generally  fixed  [Ref.  2  p. 68]  so  they  are  not  appropriate  for  the  models  de- 
veloped. A  scenario  is  set  in  Korea  for  which  realistic  data  can  be  extracted  from  avail- 
able maps  [Ref.  3j.  Consumption  data  is  provided  by  the  U.S.  Army  Ordnance,  Missile, 
and  Munitions  Center  (USAOMMCS)  [Ref.  4].  The  general  network  structure  is  ex- 
tracted from  applicable  Army  Field  Manuals  [Ref.  5  pp.2-78].  Since  specific  consump- 
tion rates  and  force  structures  tied  to  a  theater  are  classified,  our  scenario  is  deliberately 
general  and  represents  figures  which  are  realistic  but  but  not  precise. 

The  purpose  of  this  research  is  to  model  the  WADS  using  current  doctrine  to  indi- 
cate where  changes  should  be  made.  As  a  result  of  model  development,  a  tool  for  anal- 
ysis of  future  systems  is  presented  which,  with  minor  changes  to  the  source  code,  is 
complete  for  use  (assuming  availability  of  the  GAMS  software). 

C.  THE  MODEL  AND  SOLUTION  PROCEDURE 

Using  the  aforementioned  scenario,  an  integrated  approach  to  the  overall  optimiza- 
tion problem  is  described  and  preliminary  computational  results  given.  However,  results 
with  the  integrated  model  are  extremely  limited  because  it  is  too  large  for  the  available 
solver.  Consequently,  an  alternate  modelling  approach  is  described  and  analyzed  at 
length. 

A  sequential  heuristic  which  uses  a  separation  technique  solves  the  scenario  of  in- 
terest.   First,  a  binary  integer  program  locates  the  ammunition  facilities  within  a  Corps 


area.  Then,  facility  location  data  is  passed  to  a  network  flow  model  to  determine  optimal 
movement  of  ammunition  forward  from  the  Corps  rear  to  the  combat  units  in  the 
brigade  area.  The  Generalized  Algebraic  Modelling  System  (GAMS)  implements  the 
heuristic  using  a  sequential  solving  procedure. 

A  natural  approach  to  this  analysis  might  be  through  stochastic  inventory  theory. 
However,  this  approach  cannot  be  used  since  distributional  demand  data  is  not  avail- 
able. To  date,  ammunition  consumption  in  combat  has  only  been  roughly  correlated  to 
activity,  i.e.  first  day  of  defense,  first  day  of  offense,  etc.  The  Concepts  Analysis  Agency 
and  the  Combined  Arms  Center  are  presently  working  on  distribution  issues  [Ref.  6]. 
Unfortunately,  preliminary  results  are  not  available.  Therefore,  wide  use  of  combat 
simulations  provides  synthetic  data  which  is  currently  used  for  most  analysis  and  will 
be  used  in  this  thesis. 

D.     OUTLINE 

This  thesis  presents  an  integrated  and  a  sequential  approach  to  locate  ammunition 
facilities  and  determine  proper  flow.  Chapter  II  outlines  the  WADS  as  it  currently  exists 
and  describes  some  concepts  under  development  to  improve  performance  of  the  system. 
Chapter  III  develops  the  integrated  model  and  the  sequential  heuristic.  The  facility  lo- 
cation and  network  flow  portions  of  each  model  are  discussed.  Both  models  have  re- 
marks concerning  GAMS  implementation  and  highlights  of  techniques  used  to  construct 
solvable  formulations  given  sufficient  CPU  time  and  computer  memory.  Chapter  IV 
discusses  computational  experience  of  the  integrated  approach  and  the  sequential 
heuristic.  Advantages  and  disadvantages  of  the  two  procedures  are  given.  Model  be- 
havior is  outlined  which  forms  the  basis  for  conclusions  drawn  concerning  current  doc- 
trine. Chapter  V  contains  conclusions  and  recommendations  for  further  research. 
Appendices  include  the  GAMS  code  for  the  integrated  model  and  the  sequential 
heuristic. 


II.     THE  ARMY  WARTIME  AMMUNITION  DISTRIBUTION  SYSTEM 

This  chapter  outlines  the  procedures  and  structure  of  the  Wartime  Ammunition 
Distribution  System  (WADS).  The  current  system  and  future  developments  are  dis- 
cussed. Procedures  and  structure  presented  provide  the  basis  for  model  development 
discussed  in  Chapter  III. 

A.     CURRENT  AND  FUTURE  STRUCTURE 

The  WADS  [Ref.  5  :  pp. 2-78]  is  modeled  as  an  acyclic  network  directed  from  a  port, 
which  acts  as  a  source,  through  a  series  of  transshipment  nodes  (Theater  Storage  Areas, 
Corps  Storage  Areas  and  Ammunition  Supply  Points)  to  the  Ammunition  Transfer 
Points  (ATPs)  which  act  as  sinks.  Transportation  links  between  the  nodes  are  the  arcs 
of  the  network. 

Ammunition  arrives  from  the  continental  United  States  at  a  port,  within  the  TO, 
where  it  is  offloaded  from  ships  or  aircraft.  From  the  port,  ammunition  is  moved  to  the 
Theater  Storage  Area  (TSA)  or  shipped  directly  to  the  Corps  Storage  Areas  (CSAs)  and 
Ammunition  Storage  Points  (ASPs).  The  TSAs  distribute  ammunition  to  the  CSAs  and 
ASPs.  The  CSAs,  in  turn,  pass  ammunition  to  the  ASPs  and  ATPs.  ASPs  supply  am- 
munition only  to  the  ATPs.  Supply  is  based  on  a  continuous  refill  and  is  directed  from 
the  rear  to  forward  areas.  Lateral  and  forward-to-rear  movement  of  ammunition  is  not 
allowed.  Many  different  distribution  networks  are  possible,  based  on  tactical  configura- 
tions, but  a  typical  network  in  accordance  with  Army  doctrine  [Ref.  5  :  pp. 2-63]  is  shown 
in  Figure  2  on  page  6. 

In  general,  each  TSA  can  be  expected  to  support  two  or  more  CSAs.  Each  forward 
division  will  be  supported  by  at  least  one  CSA.  This  means  each  CSA  can  support  two 
or  more  ASPs.  Each  ASP,  in  turn,  could  support  one  or  more  ATPs.  For  modelling 
purposes  in  this  thesis,  each  ASP  will  support  two  ATPs  and  each  CSA  will  support  two 
ASPs  and  four  ATPs  as  indicated  in  Figure  2.   This  is  the  typical  arrangement. 

Under  current  doctrine,  stockage  of  ammunition  occurs  at  the  TSA,  the  CSA,  and 
the  ASP  on  a  major  scale,  i.e.  multiple  days  of  supply.  Stockage  at  the  ATPs  is  only 
short  term,  something  on  the  order  of  hours  versus  days.  Although  some  inconsistency 
exists  [Ref.  7  :  pp. 3-38],  the  following  table  indicates  stockage  "rules  of  thumb"  and  ap- 
proximate physical  size  [Ref.  2J: 
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Figure  2.      Current  and  Future  Wartime  Ammunition  Distribution  System 


Table  1.      SUPPLY  AND  SIZE  DATA 


FACILITY 

SUPPLY 

SIZE  (km2) 

TSA 

30  days 

20  + 

CSA 

5-7  days 

16 

ASP 

3-5  days 

9 

ATP 

3-4  hrs. 

<  <  1 

The  Army  uses  a  nebulous  concept,  "days  of  supply",  to  establish  stockage  levels. 
A  day  of  supply  varies  with  demand  and  is  defined  to  be  the  total  amount  of  ammunition 
issued  by  any  given  ATP,  and  subsequently  consumed,  during  one  day  of  combat.  For 
instance,  an  ASP  which  supports  2  ATPs  must  have,  according  to  the  previous  table, 
3-5  days  of  supply  on-hand  for  each  ATP.  A  day  of  supply  may  vary  from  ATP  to  ATP. 
ATPs  within  a  division  may  have  varying  demand  rates  and  hence  a  day  of  supply  would 
be  different  for  each.  This  means  the  2  ASPs  which  support  a  division  will  probably  have 
different  stockage  levels  since  the  days  of  supply  will  most  likely  not  be  the  same  for  each 
ATP.  Since  demand  is  random,  the  "rule  of  thumb"  for  stockage  can  fluctuate.  In  actual 
practice,  a  node  could  have  the  required  days  of  supply  on-hand  for  a  low  demand  period 
and  subsequently  violate  the  stockage  rules  in  a  matter  of  hours  should  a  period  with 
high  demand  occur.  A  day  of  supply  may  not  be  a  very  satisfactory  concept  but  it  is  the 
one  used  in  the  Army. 

By  doctrine,  each  arc  has  a  percentage  attribute  which  indicates  what  fraction  of  the 
head  node's  stockage  comes  from  the  tail  node.  These  percentages  seem  to  be  based  on 
wartime  experience  and  "best  guesses."  No  analytical  computations  appear  to  exist  that 
support  the  specified  percentages  [Ref.  6].  The  term  "bypass"  is  used  to  indicate  an  arc 
which  goes  around  and  not  through  a  transshipment  node.  Typically,  ammunition 
flowing  on  bypass  arcs  is  high  demand,  high  tonnage  items.  Low  demand,  low  tonnage 
items  usually  flow  through  each  transshipment  node  from  the  Port  to  the  ATP.  For  an 
ATP,  80%  of  ammunition  comes  directly  from  the  CSA  and  20%  from  the  ASP.  At  the 
ASP,  20%  comes  from  the  Port,  30%  from  the  TSA,  and  50%  arrives  directly  from  the 
CSA.  The  CSA  receives  50%  of  its  ammunition  from  the  Port  and  50%  from  the  TSA. 
All  of  the  TSA's  ammunition  comes  from  the  Port. 

ATPs  are  located  20-30  kilometers  (km)  from  the  forward  line  of  troops  (FLOT). 
ASPs  are  a  maximum  of  30  km  to  the  rear  of  the  ATPs  they  support  and  CSAs  are  100 


km  to  the  rear  of  the  ASPs.  This  yields  a  maximum  distance  of  130  km  from  the  CSA 
to  the  ATPs  it  supports  and  160  km  from  the  CSA  to  the  FLOT  [Ref.  5  :  pp.2-60  -  2-78]. 

Arcs  from  node  to  node  are  capacitated.  The  bulk  of  ammunition  moved  in  the  TO 
is  by  5  ton  tractors  and  22.5  ton  trailers  and  the  number  of  tractor/trailers  available  is 
finite.  Specific  upper  and  lower  bounds  for  each  arc  in  terms  of  the  ability  to  move  short 
tons  of  ammunition  by  tractor/ trailer  are  not  assigned.  The  transportation  community 
desires  maximum  flexibility  in  operations  and  prefers  not  to  dedicate  tractor/trailers  to 
individual  arcs.  By  declining  to  dedicate  assets  to  specific  arcs,  tractor/trailers  can  be 
moved  where  needed  the  most.  In  theory,  this  provides  the  best  support  [Ref.  8]. 

A  rough  transportation  plan  to  support  the  WADS  might  be  calculated  as  follows. 
For  each  time  period,  find  the  required  flows  to  meet  inventory  goals  and  demand  over 
every  arc  in  the  network.  The  determination  of  proper  flow  must  be  constrained  by 
available  transportation  assets  whose  capability  to  move  ammunition  forward  will  be 
degraded  by  maintenance  requirements  and  other  transportation  missions,  and  aided  by 
the  number  of  round  trips  possible  per  day.  The  number  of  tractor/trailers  required  for 
each  time  period  is  then  found  by  summing  all  ammunition  flows  and  dividing  by  aver- 
age haul  weight  for  each  trailer.  This  roughly  gives  the  total  number  of  tractor/trailers 
required  for  each  time  period. 

It  is  estimated  that  75-80%  of  all  cargo  moved  within  the  TO  will  be  ammunition. 
This  means  that  the  majority  of  all  transportation  assets  will  be  moving  ammunition 
forward.  Flow  forward  along  the  PORT-TSA-CSA-ASP  paths  depends  heavily  on  the 
particular  characterisitics  of  the  TO,  specifically  road,  river  and  rail  networks,  and  the 
level  of  host  nation  support,  if  any.  Therefore,  estimating  ability  to  move  ammunition 
forward  is  difficult.  On  the  other  hand,  movement  on  the  the  CSA-ASP-ATP  (Corps 
level)  paths  is  better  defined.  The  Corps  level  set  of  arcs  is  supported  almost  exclusively 
by  tactical  wheeled  vehicles;  primarily  the  5  ton  tractor  and  the  22.5  ton  trailer  whose 
average  haul  weight  is  15  tons  [Ref.  8]. 

Typically,  5  medium  truck  companies  support  a  Corps.  Each  company  is  authorized 
sixty  5  ton  tractors  and  one  hundred  and  fifty  22.5  ton  trailers.  The  transportation 
community  anticipates  an  availability  of  75%  due  to  maintenance  requirements  [Ref. 
8]. 

Normal  convoy  speeds  are  32  KiM  per  hour  for  hard  surface  roads  and  16  KM  per 
hour  for  cross  country  and  loose  surface  roads.  These  speeds  may  seem  slow,  but  for  the 
reduced  road  trafficability  typical  in  the  TO  they  are  quite  realistic. 


Each  node  within  the  WADS  has  an  ability  to  receive,  rewarehouse,  and  issue  am- 
munition. This  is  called  "lift  capacity"  (or  just  "lift")  and  is  a  function  primarily  of  ma- 
terial handling  equipment  and  personnel  assigned.  Lift  capacity  is  a  constraint  on  the 
ability  of  the  WADS  to  process  ammunition  for  inventory  or  movement  forward. 

Each  CSA  is  operated  by  one  or  more  General  Support  (GS)  Ammunition  Compa- 
nies. Each  GS  Company  has  a  current  lift  capacity  of  3696  short  tons  (STON)  per  day 
which  will  be  upgraded  to  5332  STON  per  day  in  the  future.  Equal  effort  is  usually  de- 
voted to  receipt,  rewarehousing,  and  issuing  ammunition.  This  yields  a  capacity  to  issue 
1232  STON  per  day  and  in  the  future  1777  STON.  In  high  demand  periods,  a  well 
stocked  CSA  could  use  its  lift  capacity  solely  toward  issuing  ammunition  over  some  pe- 
riod of  time  thus  increasing  the  ability  to  meet  demand.  However,  lift  is  more  properly 
thought  of  as  a  constraint  on  the  sum  of  ammunition  received,  rewarehoused  and  issued. 
By  devoting  all  lift  toward  issuing  ammunition  to  meet  a  large  demand,  receipt  of  new 
ammunition  to  replenish  stocks  or  rewarehousing  of  ammunition  on-hand  and  presently 
not  required  is  precluded. 

A  Direct  Support  (DS)  Ammunition  Company  operates  two  ASPs  which  is  the 
number  usually  allocated  to  support  a  forward  division.  The  DS  Company  has  a  current 
lift  capacity  of  2172  STON  per  day  with  a  future  upgrade  to  2732  STON  per  day.  Under 
normal  conditions,  lift  effort  is  divided  between  receipt  and  issue  of  ammunition.  Simple 
calculations  then  yield,  at  each  ASP,  a  capacity  to  issue  543  STON  per  day  and  683 
STON  per  day  in  the  future.  Under  surge  flows,  each  ASP  could  issue  1036  STON  per 
day  now  and  1366  STON  per  day  in  the  future  until  stock  exhaustion,  3-5  days  later. 
As  noted  before,  dedicating  all  lift  toward  issuing  ammunition  to  meet  demand  is  at  the 
expense  of  receiving  additional  ammunition  to  replenish  stocks.  ASPs  typically  do  not 
devote  much  effort  to  rewarehousing. 

The  number  and  capacity  of  the  ATPs  varies  with  type  of  division.  The  table  below 
indicates  different  configurations  [Ref.  5  :  p. 2-77] 

Field  testing  indicates  the  ability  of  the  ATP  to  handle  600-700  STON  per  day  for 
short  periods  of  time,  during  high  demand  periods  [Ref.  9].  Concept  papers  reviewed 
[Ref.  10]  have  projected  a  sustained  issue  of  750  STON  per  day  for  forward  ATPs  in  the 
brigade  area,  and  1450  STON  per  day  for  the  ATP  in  the  division  rear  with  appropriate 
equipment  and  personnel  upgrades. 


Table  2.      DIVISION  ATP  DATA  (F=  FORWARD,  M  =  MAIN  OR  REAR) 


Division  Type 

ATP  Density 

Capacity  per  ATP 

Armored 

Infantry 

Mechanized 

3F 

500  STON 

Heavy 

3F/1M 

350/200  STON 

Light 

3F 

250  STON 

Airborne 

3F/1M 

350/500  STON 

Air  Assault 

3F/1M 

350/500  STON 

Motorized 

3F/IM 

500/350  STON 

Port  and  TSA  capacities  will  vary  from  theater  to  theater.  Host  nation  support  units 
and  GS  Companies  would  operate  the  facilities  in  this  area  but  lift  capacities  are  not 
clear  from  available  documents  [Ref.  5  :  pp.  2-60  -  2-62]. 

Demand  for  ammunition,  seen  at  ATPs,  is  random.  Historical  data  in  the  table  be- 
low gives  average  consumption  rates  by  type  of  division  and  the  division's  current  state, 
in  STON  per  day  [Ref.  7  :  pp.  7-5  -  7-7  ]: 

Table  3.      HISTORICAL  DIVISION  CONSUMPTION 


Activity 

Type  of  Division 

Armored 

Infantry 

Mech- 
anized 

Airborne 

Air  Assault 

Defense  1st 
Day 

2432.6 
STON 

1896.3 
STON 

2156.8 
STON 

1373.4 
STON 

1825.1 
STON 

Successive 
Days 

1902.8 
STON 

1722.0 
STON 

1742.3 
STON 

1277.9 
STON 

1653.1 
STON 

Offense  1st 
Day 

1911.5 
STON 

1579.6 
STON 

1680.4 
STON 

1180.7 
STON 

1572.2 
STON 

2nd-5th 
Days 

1424.3 
STON 

1350.6 
STON 

1295.4 
STON 

1018.7 
STON 

1297.8 
STON 

Successive 
Days 

1163.4 
STON 

864.9 
STON 

1094.1 
STON 

552.2 
STON 

808.9 
STON 

Presuming  that  some  sort  of  demand  distribution  exists,  the  above  data  are  activity 
consumption  means;  variances  are  unavailable. 
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For  many  analysts,  the  above  values  are  much  too  low  and  they  have  turned  to 
simulations  for  new  consumption  rates  [Ref.  6).  The  table  which  follows  indicates  mean 
consumption  rates  for  a  heavy  division;  once  again  the  variances  are  unavailable: 

Table  4.      SIMULATED  CONSUMPTION  DATA  FOR  HEAVY  DIVISION 


Agency 

Consumption  Rate  in  STON/Day 

Logistics  Center  [Ref.  1] 

1558 

Concepts  Analvsis  and  Combined  Arms 
Center  [Ref.  1] 

2281 

USAOMMCS  [Ref.  1] 

2589 

U.S. Army-Europe  [Ref.  6] 

3470 

Others  (Upper  Bound)  [Ref.  11] 

^4545 

Significant  differences  between  historical  and  simulated  consumption  data  are  obvious 
and  have  caused  much  controversy  within  the  Army  analytic  community  [Ref.  6]. 

A  detailed  examination  of  historical  and  simulated  data  shows  artillery  is  the  major 
ammunition  consumer  comprising  65-70%  of  total  consumption  [Ref.  11]  A 
USAOMMCS  developed  list  [Ref.  12]  indicates  only  19  different  ammunition  types  ac- 
count for  91.5%  of  all  ammunition  consumed  on  the  battlefield.  Adding  one  more  am- 
munition type  to  represent  the  remaining  8.5%  of  consumption,  gives  a  useful  20 
commodity  aggregation  versus  explicitly  modeling  over  200  ammunition  types  regularly 
used.  Once  ammunition  inventory  levels  (in  STON)  are  determined  then  straight  per- 
centages based  on  the  20  major  ammunition  types  from  the  USAOMMCS  list  will  give 
a  "rough  cut"  at  individual  stockage  levels.  An  aggregated  approach  was  taken  versus  a 
multicommodity  flow  since  consumption  data  was  in  STON  not  individual  ammunition 
types,  ammunition  is  apparently  consumed  in  roughly  proportional  quantities  during  all 
time  periods,  and  the  same  material  handling  equipment  and  transportation  assets  are 
used  to  process  ammunition  for  use.  In  addition,  the  software  used  limits  the  multi- 
commodity  flow  approach  except  for  small  to  moderate  sized  models. 

B.     LOCATION 

The  essential  criterion  for  site  selection  of  an  ammunition  facility  is  close  proximity 
to  an  existing,  all  weather  road.  Distance  to  supported  units  or  storage  facilities  should 
be  minimized  within  security  constraints.  The  site  should  be  selected  to  facilitate  smooth 
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flow  of  traffic  in  and  out  once  off  the  main  road  network.  Also,  proximity  to  any  rail 
networks  is  always  welcomed  [Ref.  13]. 

Taking  advantage  of  terrain  masking  and  vegetation  is  key  to  reducing  target  sig- 
nature since  storage  facilities  and  heavy  traffic  areas  can  be  quite  obvious  if  no  effort  is 
taken  toward  concealment.  Location  away  from  other  likely  targets  is  essential  since 
ammunition  facilities  are  likely  targets  themselves. 

High  demand  units  should  be  supported  by  the  closest  facilities.  This  is  just  com- 
mon sense.  Low  demand  units  can  be  supported  by  facilities  farther  away  but  still  within 
doctrinal  limitations.  However,  Command  and  Control  requirements  may  violate  this 
common  sense  approach. 

C.     PROCEDURES 

The  Division  Ammunition  Officer  (DAO)  orders  ammunition  for  the  division  based 
on  the  Commander's  direction  and  staff  guidance.  The  order  is  determined  on  the 
knowledge  of  upcoming  operations  and  an  anticipation  of  enemy  action.  In  theory,  once 
the  DAO  places  an  order,  8-10  hours  later  ammunition  should  arrive  at  the  ATP  in  the 
quantity  and  type  ordered.  In  contrast  to  civilian  inventory  management  philosophy,  an 
ammunition  stockout  is  not  allowed. 

The  DAO  communicates  with  the  Corps  Material  Management  Center  (MMC)  who 
works  closely  with  the  Corps  Movement  Control  Center  (MCC).  The  MMC  and  MCC 
control  and  move  ammunition  from  the  Corps  rear  area  to  the  brigade  area  [Ref.  5  :  pp. 
1-27  -  1  -33 J.  The  DAO's  role  is  to  "pull"  flow  through  the  network  subject  to  MMC  and 
MCC  constraints.  The  proper  "pull"  gives  the  best  results. 

Convoys  should  arrive  at  the  ATPs  every  3-4  hours.  The  success  or  failure  of  the 
network  depends  on  the  DAO's  ability  to  keep  a  steady  flow  of  ammunition  in  the 
proper  quantity  and  type  to  meet  division  needs.  This  depends  directly  on  the  quality 
of  the  consumption  forecast  sent  to  the  MMC  and  MCC.  If  the  flow  is  too  high,  excess 
ammunition  provides  an  increased  target  signature  which  is  an  invitation  for  destruction. 
If  the  flow  is  too  low,  the  division  cannot  survive. 

The  majority  of  high  demand,  high  tonnage  items  flow  over  the  CSA-ATP  arcs 
which  are  130  km  in  length,  requiring  4-8  hours  travel  time.  Low  demand,  low  tonnage 
items  flow  over  the  ASP-ATP  arcs  which  are  30  km  long,  and  which  require  at  least  1-2 
hours  travel  time.  The  DAO's  forecast  must  be  as  accurate  as  possible.  The  DAO  can- 
not use  the  ASP  as  a  regular  source  of  ammunition,  should  his  forecast  be  continually 
in  error.  This  would  eliminate  any  emergency  source  of  ammunition  in  the  event  of  an 
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unanticipated  enemy  action.  This  all  implies  that  the  CSA-ATP  arcs  are  the  division's 
lifelines  and  the  ASP  exists  primarily  exists  for  emergency  stockage. 

D.  DISCUSSION 

The  aforementioned  network  structure  and  procedures  have  been  "tried  by  fire"  in 
World  War  II,  Korea,  and  Vietnam;  they  work.  However,  combat  under  the  Airland 
Battle  concept  will  place  an  increased  demand  on  the  system  primarily  due  to  the  in- 
creased requirements  on  the  forward  units  to  defeat  the  first  two  enemy  echelons.  Na- 
tural questions  that  have  arisen  are: 

1.  Will  the  increased  demand  for  ammunition  force  changes  in  the  current  system 
structure? 

2.  Are  the  heuristics  or  "rules  of  thumb"  from  past  wars  applicable  to  future  conflicts? 

3.  What  stockage  policies  can  be  used  to  minimize  the  size  of  ammunition  facilities 
thereby  reducing  target  signature? 

4.  What  is  the  best  placement  for  ammunition  storage  facilities  to  support  the  force? 

Logically,  one  might  question  the  need  for  the  large  number  of  storage  facilities 
(nodes)  and  numerous  bypass  (arc)  possibilities.  These  criticisms  would  be  valid  in  a 
totally  secure,  rear  area  with  numerous,  high  speed  roads  toward  the  forward  areas.  It 
is  prudent  to  assume  our  adversaries,  in  any  future  conflict,  will  disrupt  our  rear  areas 
whenever  and  wherever  possible.  Further,  road  networks  by  the  very  nature  of  war  will 
have  reduced  trafiicability  at  best  and  will  be  subject  to  interdiction.  These  two  facts 
force  a  dispersion  of  ammunition  stocks  to  reduce  target  signature  and  placement  of 
stockage  points  as  far  forward  as  possible  to  minimize  travel  time.  Redundancy  in  net- 
work flows  and  storage  points  is  necessary  to  insure  delivery  forward  should  an  arc  or 
node  be  severed  or  destroyed. 

E.  CONCEPTUAL  DIRECTIONS 

USAOMMCS,  in  an  attempt  to  best  support  the  dynamic  Airland  Battle,  is  cur- 
rently developing  the  Maneuver  Oriented  Ammunition  Distribution  System  (MOADS) 
[Ref.  10).  The  MOADS  concept  increases  the  size  of  the  CSA,  which  is  located  in  a 
relatively  secure  position,  and  reduces  the  size  of  the  ASP,  which  is  located  in  a  some- 
what insecure  position.  In  addition,  the  number  of  ASPs  serving  a  division  will  be  in- 
creased from  two  to  three.  This  will  reduce  target  signature,  increase  dispersion,  and 
provide  additional  redundancy.  Ammunition  will  be  supplied  in  Combat  Configured 
Loads  (CCLs)  for  high  demand,  high  tonnage  items.  This  will  reduce  consumer  material 
handling  requirements.  ASP  stockage  will  be  reduced  to  1-3  days  of  supply  versus  the 
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current  3-5  days  of  supply.  Ammunition  from  the  CSA  will  be  primarily  CCL  high 
tonnage,  high  demand  items.  The  ASP  will  continue  to  supply  non-CCL  low  demand, 
low  tonnage  items  and  emergency  stockage.  The  new  network  configuration  is  shown  in 
Figure  3  on  page  15. 

New  stockage  rules  for  each  ammunition  storage  facility  are  shown  below.  No  in- 
formation on  approximate  physical  size  is  available: 

Table  5.      MOADS  STOCKAGE  CRITERIA 


Facility 

Days  of  Supply 

TSA 

30 

CSA 

10 

ASP 

1-3 

ATP 

3-4  hrs. 

The  reduction  in  days  of  supply  at  the  ASP  will  reduce  the  physical  size  which  will  re- 
duce target  signature,  certainly  a  positive  benefit.  The  tradeoff  is  an  increased  target 
signature  for  the  CSA.  Since  the  CSA  is  150-160  km  to  the  rear  of  the  FLOT,  the  risk 
may  be  acceptable. 

Another  concept,  which  has  been  tested,  is  the  Palletized  Loading  System  (PLS) 
[Ref.  9].  This  idea  is  best  used  in  conjunction  with  MOADS.  Highlights  are  an  ability 
to  rapidly  upload  and  offload  flatracks  of  ammunition  pallets  thereby  reducing  material 
handling  and  personnel  requirements,  conversion  of  ATPs  to  Ammunition  Control 
Points  (ACPs)  that  direct  ammunition  convoys  to  the  unit  field  trains,  and  elimination 
of  the  TSA-ASP  bypass  arc.  The  major  advantage  is  the  rapid  upload  and  offload  of 
ammunition  which  reduces  material  handling  requirements  network-wide.  This  in  turn 
increases  lift  capacities  and  thereby  speeds  supply  to  the  consumers  forward. 

The  operational  complexity  and  manifestations  of  uncertainty  associated  with  the 
WADS  should  be  of  real  concern  [Ref.  2].  Proper  control  by  all  keys  players  can  reduce 
the  propensity  or  potential  of  the  system  to  crash.  Without  proper  control,  all  the  pro- 
posed improvements  in  the  WADS  will  be  useless. 

F.     REMARKS 

This  is  a  brief  summary  of  the  WADS.  The  current  "rules  of  thumb"  based  on  past 
conflicts  may  not  be  appropriate  for  future  wars.  Analytical  methods,  using  the  current 
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Figure  3.      MO  ADS  Wartime  Ammunition  Distribution  System 
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heuristics  as  a  starting  point,  can  provide  the  means  to  examine  the  current,  future,  and 
conceptual  systems  through  model  development,  formulation  and  solution. 
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III.     MODEL  DEVELOPMENT  AND  FORMULATION 

Two  models  of  the  WADS  are  proposed  within  this  chapter.  One  is  an  integrated 
optimization  approach  which  locates  facilities  and  determines  network  flows  concur- 
rently. The  other  is  a  separation  heuristic  which  sequentially  determines  optimal  facility 
location  and  then  flow  through  a  network  to  meet  inventory  goals  and  demand.  Both 
models  are  coded  in  GAMS  (Generalized  Algebraic  Modelling  System)  developed  by 
Meeraus  and  and  Brooke  [Ref.  14]  and  are  solved  utilizing  Marsten's  Zero-One  Opti- 
mization Methods  (ZOOM)  [Ref.  15]  or  Murtagh  and  Saunders'  Modular  In-core  Non- 
linear Optimization  System  (MINOS)  [Ref.  16]. 

A.  GAMS 

GAMS  is  a  model  generator  and  solver  interface  used  for  linear,  nonlinear,  and  in- 
teger programming.  All  work  for  model  development,  formulation,  and  solution  was 
obtained  utilizing  GAMS.  The  great  flexiblity  of  GAMS  is  its  strength.  Specifically,  the 
logical  based  "such  that"  operator,  denoted  by  the  dollar  sign  (S),  made  model  develop- 
ment, formulation,  and  subsequent  solutions  possible.  The  "such  that"  operator  is  used 
to  restrict  model  generation  to  only  those  constraints  and  variables  applicable  [Ref.  17]. 
Constraint  and  variable  reduction  made  the  difference  between  a  solvable  model  for  the 
following  scenario  and  in  some  cases  a  model  that  could  not  even  be  formulated  due  to 
size. 

B.  SCENARIO 

The  scenario  is  a  worst  case  situation  in  the  framework  of  the  WADS.  The  setting 
is  in  Korea,  where  after  a  surprise  attack  by  the  North  Korean  (NK)  forces,  we  have 
consolidated  and  are  now  vigorously  attacking  toward  the  NK  capital  of  Pyongyang. 
After  a  number  of  successful  weeks  of  attack,  we  have  stretched  our  ammunition  net- 
work to  the  limit  and  must  establish  new  ASPs  and  CSAs  forward  from  their  present 
locations.  Reconnaissance  elements  have  located  a  number  of  possible  sites  for  the  new 
ASPs  and  CSAs.  As  of  late,  NK  resistance  has  stiffened  and  an  armored  counterattack 
within  the  next  week  is  anticipated  somewhere  along  the  Corps'  FLOT.  The  Corps' 
Commander  desires  the  optimal  location  for  ammunition  facilities  utilizing  the  existing 
road  network. 
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In  addition,  the  Corps'  Commander  wonders  if  the  WADS  structure  and  current 
operating  procedures  will  be  adequate.  The  MMC  has  acquired  30  days  of  demand  data 
from  another  Corps'  past  consumption  data  under  similar  circumstances.  The  Corps' 
Commander  would  like  an  analysis  of  the  Corps'  WADS  based  on  this  data. 

The  Corps  is  composed  of  three  Divisions  abreast,  in  contact,  and  the  Corps  rear 
units  which  support  the  forward  Divisions.  Each  Division  will  have  4  ATPs  which  gives 
the  Corps  12  ATPs  total.  Each  Division  will  be  supported  by  2  ASPs  and  1  CSA.  This 
means  that  the  total  Corps'  requirement  for  ASPs  is  6  and  for  CSAs  is  3.  A  total  of  12 
possible  ASP  sites  have  been  declared  suitable  and  4  tentative  CSA  locations  selected. 
The  existing  road  infrastructure  provides  the  network  and  distances  between  facilities  is 
road  distance  unless  stated  otherwise.  The  idea  is  gradually  to  close  down  the  old  ASPs 
and  CSAs  to  the  rear  and  open  new  sites  forward.  A  graphical  representation  of  the 
tactical  situation  is  shown  in  Figure  4  on  page  19. 

The  Corps'  three  Divisions  will  be  composed  of  two  Mechanized  Divisions  posi- 
tioned on  the  flanks  and  one  Armored  Division  in  the  middle.  Cynics  will  quickly  point 
out  that  mechanized/armored  forces  are  not  appropriate  for  the  Korean  theater.  How- 
ever, keeping  in  sight  the  goals  of  a  generic  Corps  not  tied  to  a  specific  theater  and  an- 
alyzing the  structural  aspects  of  the  WADS,  consistency  is  maintained. 

C.    THE  INTEGRATED  MODEL  DEVELOPMENT 

The  integrated  approach  solves  the  ammunition  facility  location  and  ammunition 
flow  problems  concurrently.  This  is  accomplished  by  using  a  generalized  form  of  the 
capacitated  facility  location  model  [Ref.  18  :  pp.  195-197]  which  is  a  mixed  integer  pro- 
gram (MIP).  The  capacitated  facility  location  model  places  upper  bounds  on  the  capac- 
itates of  facilities.  In  the  problem  at  hand,  the  model  places  upper  bounds  on  the  ability 
of  the  ammunition  storage  facilities  to  issue,  rewarehouse  and  receive  ammunition.  Up- 
per and  lower  bounds  are  established  for  inventory.  The  doctrinal  support  relationships 
between  ammunition  storage  facilities  outlined  in  Chapter  II  and  an  underlying  distrib- 
ution network  with  side  constraints  (based  on  "rules  of  thumb"  developed  in  past  con- 
flicts) make  this  problem  difficult  to  solve.  Binary  variables  determine  ammunition 
facility  location  and  support  relationships  between  storage  facilities.  Continuous  vari- 
ables are  ammunition  flow  from  node  to  node  and  inventory  from  time  period  to  time 
period.  Binary  and  continuous  variables  interact  in  the  constraint  matrix  where  condi- 
tional constraints  are  used.  If  solved,  the  MIP  gives  an  optimal  solution  with  respect  to 
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Figure  4.      Scenario 
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facility  location  and  network  flow.  This  would  be  the  most  natural  way  to  solve  the 
scenario  if  computationally  feasible. 

D.     THE  INTEGRATED  FORMULATION 

The  entire  scenario  is  modeled  as  a  generalized  capacitated  plant  (ammunition  stor- 
age facility)  location  problem  which  is  a  mixed  integer  program  (MIP).  The  following 
formulation  will  open  (1)  or  close  (0)  facilities  using  binary  variables.  Continuous  vari- 
ables are  ammunition  flows  between  facilities. 

1.  Indices 

/=  1,2,...,/  (=  12)  ATPs  (fixed) 
j  —  1,2,...,/  (=  12)  ASPs  (proposed  sites) 
k  =  1,2,...,AT  (  =  4)  CSAs  (proposed  sites) 
t  =  l,2,...,r ( =  30)  time  periods 

2.  Data 

dlt   demand  in  STON/ 1000  at  ATP  i  in  time  period  t 

/(/   distance  by  road  from  ATP  i  to  ASP  j  in  km 

ljk  distance  by  road  from  ASP  j  to  CSA  k  in  km 

ltk  distance  by  road  from  ATP  i  to  CSA  k  in  km 

pv  penalty  cost  for  road  quality  from  ATP  i  to  ASP  j 

pjk   penalty  cost  for  road  quality  from  ASP  j  to  CSA  k 

pik   penalty  cost  for  road  quality  from  ATP  i  to  CSA  k 

N    number  of  ASPs  per  Corps  ( =  6) 

M    number  of  CSAs  per  Corps  (  =  3) 

clJt  cost  to  move  one  STON  of  ammunition  from  ASP  j  to  ATP  i  in  time  t 

cJkt   cost  to  move  one  STON  of  ammunition  from  CSA  k  to  ASP  j  in  time  t 

crt,  cost  to  move  one  STON  of  ammunition  from  CSA  k  to  ATP  i  in  time  t 

hJt  cost  to  hold  one  STON  of  ammunition  at  ASP  j  from  time  t-I  to  t 

ML"'   maximum  lift  capacity  in  STON/ day  of  the  CSA  ( =  10664) 

ML"  maximum  lift  capacity  in  STON/day  of  the  ASP  (=  2732) 

P™*   percentage  of  demand  provided  by  the  ASP  ( =  0.2) 

/**•  percentage  of  demand  provided  by  the  CSA  ( =  0.8) 

/•»  maximum  number  of  days  of  supply  at  the  ASP  ( =  5) 

/'••   minimum  number  of  days  of  supply  at  the  ASP  (=1) 

5mln    minimum  distance  from  ATP  i  to  the  FLOT  in  km  (  =  20) 

Sj  distance  from  ASP  j  to  the  FLOT  in  km 
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sk   distance  from  CSA  k  to  the  FLOT  in  km 

u   adjustable  constant  for  the  tactical  situation  (  =  1) 

a     scaling  factor  for  the  objective  function  ( =  .01) 

v     number  of  round  trips  per  day  from  CSA  k  to  ATP  i  ( =  3) 

r     number  of  round  trips  per  day  from  CSA  k  to  ASP  j  (  =  4) 

T    number  of  tractors  authorized  per  Corps  (=  300) 

h     average  haul  weight  per  trailer  in  STON/1000  (=  .015) 

a     tractor/trailer  availability  on  an  given  day  (  =  0.8) 

pammo    percentage  of  total  tractor/trailers  hauling  ammunition  (  =  0.8) 

w     work  level  at  CSA  issuing  ammunition  ( =  .333) 

C,   cardinality  of  the  index  set  t 

3.  Decision  Variables 

x,j    1  if  ATP  i  is  supported  by  ASP  j,  0  otherwise 
aJk    1  if  ASP  j  is  supported  by  CSA  k,  0  otherwise 
zlk    1  if  ATP  i  is  supported  by  CSA  k,  0  otherwise 
y}    1  if  ASP  is  located  at  site  j,  0  otherwise 
bk    1  if  CSA  is  located  at  site  k,  0  otherwise 
f,Jt  flow  from  ASP  j  to  ATP  i  in  time  period  t 
fJkt   flow  from  CSA  k  to  ASP  j  in  time  period  t 
flkl  flow  from  CSA  k  to  ATP  i  in  time  period  t 
lj,    inventory  at  ASP  j  at  the  end  of  time  period  t 

4.  Formulation  (GAMS  equation  designations  in  brackets) 

/       /  J       K  IK 

min  2^  2^a  hj  Pij  Xy  +  Y  \<x  lJk  pJk  aJk  +  Y  V a  llk  plk  zlk  + 


J         K         l  I         K.         1  J         1 


[  OBJFCN  ] 
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Subject  To: 
2/fr-l     V/  (2.2) 


J 


J- 


[ONESITE(/)] 

Yjr„-ty     V/  (2.3) 

[  SERVICE  0")] 
fy^ML^Xy     Vi,j,t  (2-4) 

[  BOUND  1  (/,  y,  r)] 

^<^y      V/,y  (2.5) 

[  VARUPBD1  (i,j)] 

K 

YaJk=yj     V/  (2.6) 

[HANNA(y)] 


I 


[  SUPPORT  (  k  )] 
Jj*  £AfI"%    vy,*,/  (2-8) 

[BOUND2  0',  A,  f)] 
*,*<:**     V/,£  (2.9) 
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[  VARLPBD2  (j,  k  )] 


z,b  =  1     Vi 


-Ik 


(2.10) 


[  SINGLE  (/)] 

/ 


;=i 


[HELP  (A)] 
flkt£MLcsazlk     VU,f 

[  BOUND3  ( /,  A,  / )] 
2ik<>bk     V/,  A 

[  VARUPBD3  ( /,  A  )] 


7=1 


[  LIMIT] 


(2.11) 


(2-12) 


(2.13) 


(2-14) 


A  =  M 


(2-15) 


[  CEILING  ] 


/<=i  (=i 


[BALASPU  r)] 
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(2.16) 


J  A 


7=1  F=l 


(2.17) 


[  BALATP(/,  r)] 


/=1  (=1 


(2.18) 


[CAPCSA(£,  0] 


\fm+YJf<J^MLaspyj  v*' 


<=\  i=\ 


(2.19) 


[CAPASP(y,  0] 


J       T  IT 


(2.20) 


[STABLE  (A)] 


J       K  IK 

v^  v^  fjkl   ,  V  V  ^ 


LL     a    +  Z.Z,     a    ^ThaPammo     V/ 


(2.21) 


[  TRANS  (f)] 


I'™  pk  r+1  *<,  ^  lJt  Z  IUP  Ydit+}  Xy       V/,  t  =  l,2,...,r-  1 


(2.22) 


[  LOWER  and  UPPER  (j,  t )] 


/     a:     t 


i     r 


EZZa.-'-ZI* 


i=l   £=1   f=l 


i'=l  ;=l 


(2.23) 


[ LONGRUN ] 
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fjktJijffiktJjt^O  (2-24) 

xij>ajk'zik>yj>bke  I0.1)  (2-25) 

The  GAMS  code  for  the  MIP  follows  in  Appendix  A. 

Decision  variables,  x,;  and  aJk,  establish  a  support  relationship  so  ammunition 
may  flow  from  CSA  k  through  ASPy  to  ATP  /.  Decision  variable  zlk  allows  ammunition 
to  flow  on  the  bypass  arc  from  CSA  k  to  ATP  i. 

The  objective  function,  equation  2.1,  is  composed  of  two  parts,  one  for  locating 
the  ammunition  facilities  and  the  other  for  flow  of  ammunition. 

Lengths  /,,,  lJk,  and  llk,  and  penalty  costs  ptJ,  pjk  ,  and  plk,  form  the  basis  for  one 
portion  of  the  objective  function.  Essentially,  the  penalty  costs  make  road  distance 
longer  for  use  of  substandard  roads  (arcs).  In  particular,  a  cost  of  1.00  is  given  to  two 
lane  roads,  2.00  for  one  lane  roads,  and  4.00  for  trails.  For  composite  roads,  a  linear 
combination  is  used. 

To  fix  a  cost  for  direct  distance  from  the  FLOT  to  the  ATPs,  proximity  costs 

5  5 

were  established  utilizing  a  convex  cost  function  ( — t2-  )"  and  ( — ~ )"  where  smin  is  mini- 
mum  distance  from  the  ATP  to  the  flot.  The  values  s,  and  sk  are  the  straight  line  dis- 
tances from  the  ASP  and  CSA  to  FLOT.  Multiplying  by  s,  scales  the  convex  function  in 
terms  of  the  ATP  and  u  is  an  adjustable  parameter  which  may  vary  with  the  tactical 
situation.  As  the  distance  to  the  FLOT  decreases,  the  cost  increases  which  will  discour- 
age placement  of  an  ASP  or  CSA  too  far  forward  where  it  might  be  destroyed.  This 
portion  of  the  objective  function  forms  the  fixed  costs  for  placing  an  ASP  at  site  7  and 
placing  a  CSA  at  site  k. 

The  trade-off  between  length  modified  by  a  road  penalty  cost  and  the  convex 
function  based  proximity  cost  "drive"  the  location  portion  of  the  objective  function.  The 
modified  length  portion  minimizes  distance  from  facility  to  facility  and  has  a  tendency 
to  "pull"  the  ASPs  and  CSAs  toward  the  FLOT.  The  fixed  proximity  costs  "push"  the 
ASPs  and  CSAs  away  from  the  FLOT.  The  adjustable  scalar  a  regulates  the  amount  of 
"push"  and  "pull." 

The  flow  portion  of  the  objective  function  is  straightforward.  Costs  per  unit 
c,jr,  cJkn  and  cikt  are  charged  for  moving  ammunition  while  hJt  is  charged  for  holding  each 
unit  of  ammunition.  These  costs  are  subjective  and  are  determined  by  the  user.  For  in- 
stance, the  costs  might  be  thought  of  a  risk,  and  therefore  would  change  with  the  bat- 
tlefield situation  and  time. 
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Many  of  these  constraints  are  conditioned  on  whether  a  ammunition  storage 
facility  is  open  /  closed  or  support  relationship  between  storage  facilities  is  established  / 
not  established.  If  an  ammunition  facility  is  open  or  support  relationship  established, 
the  following  constraints  are  "turned  on":  2.3,  2.4,  2.5,  2.6,  2.7,  2.8,  2.9,  2.11,  2.12,  2.13, 
2.18,  2.19,  2.20  and  2.22.  Otherwise,  when  the  ammunition  storage  facility  is  closed  or 
the  support  relationship  not  established,  the  aforementioned  constraints  left  hand  side 
and  right  hand  side  are  set  to  zero. 

The  description  of  each  constraint,  below,  is  the  result  of  doctrine  or  required 
to  formulate  a  solvable  model.  Upper  bounds  on  lift  capacity  for  the  CSAs  and  ASPs 
as  well  as  upper  and  lower  bounds  for  inventory  are  outlined  in  Chapter  II.  Constraint 
2.2  specifies  only  one  support  path  from  ATP  /  to  ASP  j  may  exist  for  each  ATP  /. 
Constraint  2.3  forces  ASPy  to  support  two  ATPs,  if  opened.  Constraint  2.6  establishes 
a  unique  support  relationship  between  CSA  k  and  ASP  j,  if  opened.  Constraint  2.7 
mandates  CSA  k,  if  opened,  to  support  two  ASPs.  Constraint  2.10  specifies  that  only 
one  CSA  k  may  support  a  given  ATP  i.  Constraint  2.11  requires  a  CSA,  if  opened,  to 
support  four  ATPs.  Constraints  2.14  and  2.15  set  the  Corps  authorization  for  ASPs  and 
CSAs  respectively.  Constraints  2.16  and  2.17  are  flow  balance  equations  for  the  ASP  and 
ATP:  ammunition  flow  in  must  equal  flow  out.  Constraint  2.18  requires  all  flow  out  of 
CSA  k  to  ASPy  or  ATP  /,  not  to  exceed  its  lift  capacity.  Constraint  2.19  insures  that  the 
amount  of  ammunition  received  and  issued  by  the  ASP  does  not  exceed  the  established 
lift  capacity.  Constraint  2.20  sets  the  upper  and  lower  bounds  for  inventory  at  any  ASP, 
if  opened.  Constraint  2.21  establishes  a  long  run  upper  bound  of  33%  on  the  amount 
of  lift  capacity  that  a  CSA  may  devoted  to  issuing  ammunition.  This  is  a  realistic  re- 
quirement since  a  CSA  must  devote  lift  effort  not  only  to  issuing  ammunition  but  re- 
ceiving new  ammunition  and  rewarehousing  ammunition  on-hand.  Corps  transportation 
assets  form  the  basis  for  constraint  2.22.  Essentially,  flows  out  of  all  CSAs  to  each  ASP 
and  ATP,  divided  by  round  trips  possible  per  day,  must  be  less  than  available  transpor- 
tation (degraded  by  maintenance  and  mission  factors)  times  average  haul  weight  per 
trailer.  By  doctrine,  Corps  transportation  assets  move  ammunition  from  the  CSAs  to  the 
ASPs  and  ATPs.  Division  transportation  assets  move  ammunition  from  the  ASPs  to  the 
ATPs.  Finally,  constraint  2.23  supports  doctrine,  which  requires  the  majority  of  am- 
munition supplied  to  an  ATP  to  come  from  the  CSA  on  bypass  arcs. 

Constraints  2.4,  2.5,  2.8,  2.9,  2.12,  and  2.13  are  called  variable  upper  bounds 
(VUBs)  by  Schrage  [Ref.  18  :  pp.  193-195].  The  VUBs  create  a  much  tighter  formulation 
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of  the  MIP,  often  resulting  in  natural  integer  solutions  [Ref.  19  :  pp.6 1-68].  Constraints 
2.5,  2.9,  and  2.13,  in  the  context  of  this  scenario,  state  that  if  AS?  j  or  CSA  k  are  not 
open  then  no  support  paths  xu,  a]h,  zlk  are  possible.  Constraint  2.4  allows  a  flow  of  am- 
munition from  ASP  j  to  ATP  i,  if  opened,  and  establishes  an  upper  bound  on  flow  given 
by  the  lift  capacity  of  the  ASP  during  time  period  t.  Constraint  2.8  allows  flow  from 
CSA  k  to  ASP  j  if  a  support  relationship  is  established  and  once  again  sets  an  upper 
bound  on  flow.  The  lift  capacity  of  the  ASP  is  used  to  bound  flow  since  ammunition 
flowing  from  CSA  k  to  AS?  j  during  time  period  t  must  not  exceed  the  capability  of  the 
ASP  to  receive  that  flow.  Constraint  2.12  sets  an  upper  limit  of  flow  on  the  bypass  arcs 
if  a  support  relationship  exists  between  CSA  k  and  ATP  i.  Lift  capacity  of  the  CSA  sets 
the  upper  bound  on  flow  to  ATP  /  in  time  period  /.  This  simple,  yet  elegant  concept  of 
the  VL'Bs  is  extremely  powerful  as  it  often  avoids  the  requirement  for  branch  and  bound, 
cuts,  or  a  heuristic  to  arrive  at  an  integer  solution  [Ref.  18  :  pp. 193-195].  The  drawback 
is  the  large  number  of  constraints  generated  which  increases  with  index  size,  i.e  x1Jk  <>}>; 
has  I  x  J  x  K  constraints.  Since  optimal  integer  solutions  for  the  MIP  are  difficult  to 
obtain  without  the  VUBs,  the  computational  burden  they  impose  must  be  borne  [Ref. 
19  :  pp.61-68]. 

The  power  of  the  GAMS  modelling  language  is  easily  seen  as  thousands  of 
constraints  are  generated  with  only  42  equation  generation  commands!  The  amount  of 
effort  required  to  formulate  this  model  in  say  LINDO  [Ref.  18]  would  be  significantly 
greater. 

The  integrated  formulation,  although  mathematically  correct,  is  not  able  to  ar- 
rive at  feasible  solutions  (in  a  reasonable  period  of  time)  for  the  problem  at  hand. 
Therefore,  another  approach  is  required. 

E.     THE  SEQUENTIAL  HEURISTIC  DEVELOPMENT 

Because  full-scale  solutions  to  the  MIP  proved  to  be  impractical,  a  sequential 
heuristic  was  developed  which  separates  the  problem  into  one  of  first  locating  ammuni- 
tion facilities  and  then  establishing  valid  flows  between  these  facilities.  The  heuristic 
takes  advantage  of  a  sequential  solving  procedure  available  to  GAMS  and  also  follows 
the  principal  of  only  entering  data  once. 

Part  one  of  the  heuristic  is  a  facility  location  problem.  This  problem  is  modeled  as 
a  binary  integer  program  (BIP)  known  in  the  literature  as  the  uncapacitated  plant  lo- 
cation model  [Ref.  18  :  pp. 193-195]  or  the  m-median  problem  [Ref.  19  :  pp. 58-60].   This 
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BIP  can  be  extended  to  optimally  locate  ammunition  storage  facilities  in  accordance  with 
criteria  outlined  in  Chapter  II. 

Part  two  of  the  scenario  can  be  modeled  as  a  minimum  cost  network  flow  problem, 
specifically  a  production-transportation-inventory  (PTI)  model  [Ref.  20].  In  this  sce- 
nario production  will  be  called  supply  and  a  supply  of  ammunition  will  flow  through  the 
network.  In  its  simplest  form,  the  PTI  model  is  a  pure  network.  However,  due  to  doc- 
trinal restrictions,  side  constraints  exist  that  cannot  be  eliminated.  This  is  key  since  with 
a  pure  network,  specialized  solvers,  such  as  GNET  [Ref.  21],  could  be  used  to  determine 
optimal  flow.  Unfortunately  the  network's  side  constraints  preclude  the  use  of  GNET 
and  GAMS'  LP  solver  is  utilized  instead.  Solution  times  will  be  slower  than  a  specialized 
solver  but  still  acceptable. 

Solving  the  BIP  and  PTI  models  in  succession  through  a  GAMS  sequential  solving 
technique  allows  the  development  of  a  heuristic  that  separates  the  problem  into  2  small 
submodels.  The  BIP  is  solved  and  facility  solution  decisions  passed  to  the  PTI  model  to 
determine  flows. 

F.     THE  SEQUENTIAL  HEURISTIC  FORMULATION 

The  sequential  formulation  differs  not  only  in  the  separation  technique  used  but  in 
the  relationship  between  binary  (location)  and  continuous  (flow)  variables.  In  the  MIP, 
location  and  flows  are  related  together  in  the  constraint  matrix.  In  the  sequential 
heuristic,  this  relationship  is  reflected  only  in  the  objective  function  of  the  BIP  by  a  de- 
mand term  (which  must  be  satisfied  by  flow).  The  BIP  contains  only  binary  variables. 
The  PTI  uses  continuous  variables. 

Briefly,  the  facility  location  problem  is  solved  to  determine  those  ASPs  and  CSAs 
opened  and  support  relationships.  Next,  the  binary  data  is  molded  into  "such  that"  op- 
erators which  eliminate  network  constraints  that  are  not  applicable.  Finally,  a  small 
network  flow  model  is  solved  resulting  in  flows  necessary  to  meet  inventory  goals  at  the 
ASPs  and  demand  at  the  ATPs. 

The  sequential  heuristic  is  formulated  as  follows: 

1.  Indices  (same  as  MIP) 

2.  Data  (same  as  MIP) 

3.  Binary  Decision  Variables  for  the  BIP 

xlJk    1  if  ATP  i  is  supported  by  ASP  j  which  is  supported  by  CSA  k,  0  otherwise 
zlk    1  if  ATP  i  is  supported  by  CSA  k,  0  otherwise 
ys    1  if  ASP  is  located  at  site  j,  0  otherwise 
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bk    1  if  CSA  is  located  at  site  k,  0  otherwise 

4.  Continuous  Decision  Variables  for  the  PTI 

fiJt   flow  from  ASP  j  to  ATP  i  in  time  period  t 
fjkt   flow  from  CSA  k  to  ASP  j  in  time  period  t 
fikt   flow  from  CSA  k  to  ATP  i  in  time  period  t 
Ij,   inventory  at  ASP  j  at  the  end  of  time  period  t 

5.  Formulation  (GAMS  equation  designations  in  brackets) 

The  first  module  is  the  uncapacitated  facility  location  problem  which  is  a  B1P. 

i      J      K      T 


min  2_j  A  A  Z£*  liJ  PiJ +  a  hk  Pjk\  Pasp  dit  Xljk  + 


i=\  j=\  k=\  r=l 


/       K       T 


« ht  p*  p"a  <*«  * + X(  ^T  )U» + Z(  ^ir )U  bk  (126) 


1=1  j=i  ;=„ 


[  OBJFCN  ] 
Subject  To: 


J       K 


£2/«k-l     Vi  (2.27) 


y=i  k=\ 


I       K 


(=1  k^l 


I      J 


[ONESITE(/)] 

K 

Y^^=2yj  v  (2-28) 

[SERVICER)] 

j 
XX^  =  4**     VA  (2.29) 

/=1  y=l 

[  SUPPORT  (/c)] 
Jty*^     V/,y,A  (2.30) 
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[  VARUPBD1  (i,j,  k)-] 
xyk^bk     VlJ,k  (2.31) 

[  VARUPBD2  ( /,  j,  k  )] 

K 

Y^=l     V/  (2.32) 


k=\ 


[  SINGLE  (/)] 

Yz^Abj;     V*  (2.33) 


[HELP(*V] 
zik<,bk     V/,A  (2.34) 

[  VARUPBD3  ( /',  rt  )] 
y 

7=1 

[  LIMIT] 

K 

Vbk  =  A/  (2.36) 

[CEILING] 
*(#>*&..#  A  e  {0,1}  (2.37) 

The  binary  solution  data  is  now  passed  to  form  "such  that"  operators  for  the 
second  module.  The  "such  that"  operator  (S)  forms  the  following  sets: 
J'  =  ASPs  which  are  open. 
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K  =  CSAs  which  are  open. 
L{  I,J')  =  ASPj  supports  ATP  i 
L{  J',K')=  CSA  k  supports  ASPj 
L{  /,  K' )  =  CSA  k  supports  ATP  / 

The  network  flow  model  (PTI)  is  now  formulated  with  doctrinal  side  constraints  over  the 
sets  defined  above: 

T  T 

™n       Z       Xc</'^+        Z        Xw*r  + 

(i,y')eL</,/')  '=1  (j.k)eL(J ',  K' )  f=l 

[  NOBJFCN  ] 

r  r 

X    Zc^+ZZvyv  (2-38) 

(l,k)eL(l,  K')  f=l  ye/'  f=l 

Subject  To: 
Z        -&+       Z       4-^-1  +  ^  =  0    V/e7',r  (2.39) 

U,k)*L(J',K')  (i,j)<=L(l,J') 

[BALASPU  0] 
"        Z       -fr-         Z        Ar  — ^     W'<  (2-4°) 

(/J)eL(/,/')  (l,k)eL(I,K') 

[BALATP(/,  f)] 

Z         fjkt+         Z        fikt^MLcsa     VkeK',t  (2.41) 

=L(/',  ff')  (',*) 

[CAPCSA(A,  f)] 


Z         ■&  +        Z       fv^MLaSP    y/e/'.r  (2.42) 

(J.k)mUJ',  K')  (IJ)eL(I,  J') 
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U,k)eL(J',  K')  (l,k)eL(I,  K') 


[CAPASPC/,  /)] 

T  T 

Z    Z^<+    Z    Z^"c<A/z/ja  ykeK>  (2-43^ 


U,k)eL(J',K')  t=l  (i,k)eL(I,K')  i=\ 


[  STABLE  (*)] 


I 


—  4-         V  Ik.^ThaPammo     W  (2.44) 


J  A 

(J,k)eL(J',K')       V  (l,k)eL(l,K')      r 


[  TRANS  (r)] 

r  /     r 

Z   Z^=p"ZZ^  <2-45) 

(t,kMl,K')  i=l  (=1   t=\ 

[ LONGRUN  ] 
/^XX'+i  <  /;,^  /"'/>//,+,     V/e  7 ',  r  =  1,2,...,7M  (2.46) 

[  Defined  in  Variable  Declaration  ] 
fjkofijtJikt.IjtZO  (2.47) 

The  GAMS  code  for  the  sequential  heuristic  follows  in  Appendix  B. 

The  sequential  heuristic  formulation  is  quite  similar  to  the  MI  P.  In  fact,  many 
of  the  same  constraints  are  used.  The  major  difference  is  the  decision  variable  xIJk  instead 
of  x,j  and  aJk.  The  variable  xiJk  establishes  a  single  support  path  from  the  CSA  k  through 
ASP  j  to  ATP  i  versus  two  paths  with  decision  variables  xtJ  and  aJk. 

The  objective  function's  structure  for  the  first  module  is  different  from  the  lo- 
cation portion  of  the  MIP.  Here  a  trade- off  between  demand  and  length  modified  by  a 
road  penalty  is  created.  Those  ATPs  with  higher  demands  aggregated  over  T  time  peri- 
ods will  be  supported  by  the  ASPs  and  CSAs  which  are  the  closest.  Those  ATPs  with  less 
demand  over  T  time  periods  will  be  supported  by  ammunition  facilities  farther  away. 
Since  the  majority  of  demand  is  supplied  over  the  CSA-ATP  arcs,  demand  is  weighted 
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by  the  long  term  bypass  goal  (80%)  for  the  CSA  positioning.  ASP  location  is  determined 
by  the  remaining  weighted  demand.  Incorporating  demand  in  the  objective  function  is 
necessary  for  the  sequential  heuristic  since  flow,  inventory  and  location  variables  do  not 
interact  in  the  constraint  matrix. The  decomposition  technique  separates  location  from 
flow  and  inventory. 

Constraints  2.27  through  2.37  generate  the  constraint  matrix  for  the  first  mod- 
ule, an  uncapacitated  facility  location  problem,  a  B1P.  Constraint  2.27  specifies  only  one 
ijk  support  path  may  be  established  to  ATP  /.  Constraint  2.28  forces  ASPy  to  support 
two  ATPs  /.  Constraint  2.29  allows  for  two  ASPs  J  to  be  directly  supported  by  CSA  k 
and  four  ATPs  /  indirectly.  Constraint  2.32  establishes  only  one  ik  bypass  from  CSA  k 
to  ATP  /.  Constraint  2.33  mandates  CSA  k  to  support  four  ATPs  /  directly.  Constraints 
2.35  and  2.36  set  the  Corps  authorization  for  ASPs  and  CSAs  respectively.  Constraints 
2.30,  2.31  and  2.34  are  VUBS. 

After  solving  the  BIP,  the  binary  data  which  establishes  support  relationships 
and  those  ASPs  and  CSAs  opened  is  passed  through  an  error  check.  If  no  errors  are 
detected,  then  "such  that"  operators  are  coded  to  generate  the  exact  number  of  con- 
straints necessary.  For  the  problem  at  hand,  over  70,000  are  possible.  Logical  variable 
and  constraint  elimination  resulted  in  over  a  98%  reduction  to  1332  constraints,  well 
within  the  scope  that  ZOOM  or  MINOS  can  solve  [Ref.  14  :  Ch.18  pp.3-4]. 

The  objective  function  for  the  second  module  is  exactly  the  same  as  the  flow 
portion  of  the  MIP,  i.e.,  cost  to  move  and  cost  to  hold.  Constraints  are  similar  except 
that  only  constraints  corresponding  to  the  open  facilities  in  the  BIP  solution  are  allowed. 
Constraints  2.39  and  2.40  are  flow  balance  constraints  for  the  ASPs  and  ATPs  respec- 
tively. Constraints  2.41,  2.42,  2.43,  2.44  and  2.45  disrupt  the  pure  network  flow.  Con- 
straints 2.41  and  2.42  capacitate  the  the  CSAs  and  ASPs.  Constraint  2.43  sets  a  long 
term  goal  on  percentage  of  lift  capacity  devoted  to  issuing  ammunition.  Constraint  2.44 
insures  flow  does  not  exceed  transportation  capabilities  and  constraint  2.45  establishes 
doctrinal  bypass  requirements.  Constraint  2.46  provides  an  upper  and  lower  bound  for 
the  inventory  at  ASP  j.  By  doctrine,  a  required  number  of  days  of  supply  are  to  be  on- 
hand  for  each  ATP  an  ASP  supports. 

G.    TRANSPORTATION 

Required  transportation  assets  to  support  ammunition  flow  are  determined  at  the 
end  of  each  GAMS  formulation.  Ammunition  flows  are  constrained  by  available  trans- 
portation assets  which  are  degraded  and  aided  by  the  factors  discussed  in  Chapter  II. 
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Constraint  2.21  in  the  MIP  and  constraint  2.44  in  the  sequential  heuristic  restrict  am- 
munition flows  from  CSA  k  to  ASP  j  and  ATP  /  in  accordance  with  number  of 
tractor/trailers  available,  average  haul  weight  per  22.5  ton  trailer,  availability  due  to 
maintenance  requirements  and  percentage  of  Corps  transportation  assets  dedicated  to 
moving  ammunition  forward.  Ammunition  flows  between  storage  facilities  are  divided 
by  number  of  round  trips  possible  per  day  which  levels  out  flow  over  a  24  hour  period. 
The  number  of  5  ton  tractors  and  22.5  ton  trailers  needed  is  calculated  by  dividing  the 
restricted  flows  with  an  average  haul  weight  of  15  STON.  A  rough  transportation  plan 
is  then  obtained  which  gives  the  required  number  of  tractor/trailers  for  each  set  of  arcs 
to  move  ammunition  forward  to  meet  inventory  goals  and  demand.  This  estimate  is  no 
doubt  optimistic  since  it  assumes  that  tractor/ trailers  are  an  infinitely  divisible  resource. 
An  argument  could  be  made  to  include  a  relocation  factor  in  the  right-hand-side  of 
constraints  2.21  and  2.44  to  account  for  the  optimism  in  calculations. 
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IV.     COMPUTATIONAL  EXPERIENCE  AND  SOLUTIONS 

This  chapter  presents  computational  results  and  scenario  solutions.  Implementation 
of  the  .VHP  formulation  and  results  for  a  3  time  period  problem  are  discussed.  Results 
for  the  sequential  heuristic  are  then  reported  and  a  comparison  of  both  the  integrated 
and  sequential  approaches  is  given.  Finally,  the  chapter  is  closed  with  comments  con- 
cerning the  utilization  of  GAMS  for  modelling  the  WADS. 

A.     THE  MIP 

The  MIP  for  the  given  30  time  period  scenario  could  not  be  formulated  on  the 
computer  due  to  the  size  of  the  problem.  The  only  option  before  abandoning  the  MIP 
in  favor  of  the  sequential  heuristic  was  to  reduce  the  number  of  time  periods  until  the 
problem  was  small  enough.  A  reduction  from  30  time  periods  to  3  time  periods  was  re- 
quired before  a  formulation  and  solution  could  be  obtained. 

Using  an  IBM  3033  AP  for  a  3  time  period  formulation  of  the  MIP,  ZOOM  requires 
86.6  CPU  seconds  to  arrive  at  an  optimal  solution.  The  LP  relaxation  of  the  MIP  gives 
an  integer  solution. 

For  the  sake  of  comparison,  the  sequential  heuristic  was  formulated  for  3  time  pe- 
riods. Using  ZOOM  for  the  first  module  and  MINOS  for  the  second  module  resulted  in 
solution  times  of  16.9  CPU  seconds  for  the  ammunition  facility  location  model  and  0.5 
CPU  seconds  for  the  network  flow  model.  This  yields  a  total  of  17.4  CPU  seconds, 
roughly  5  times  faster  than  the  MIP! 

The  solution  values  for  the  MIP  and  sequential  heuristic  differ  in  ASPs  opened.  The 
MIP  gives  the  best  solutions  if  it  can  be  solved.  Solution  times  and  the  size  of  problems 
that  can  be  solved  certainly  favor  the  sequential  heuristic.  However,  an  examination  of 
solution  quality  is  necessary  to  fully  justify  general  use  of  sequential  heuristic  for  analy- 
sis. 

By  fixing  the  ASP  and  CSA  open/close  values  at  solution  levels  given  by  the  se- 
quential heuristic  and  forcing  those  values  into  the  MIP,  the  quality  of  solutions  can  be 
measured  by  the  change  in  the  objective  function  value.  A  small  change  in  objective 
function  value  would  indicate  that  the  sequential  heuristic  solutions  are  "good"'  and 
would  justify  use  in  the  following  analysis. 

The  objective  value  of  the  3  time  period  MIP  from  a  "cold  start"  is  81.75.  The  LP 
relaxation  gives  an  integer  solution.  Using  the  ASP  and  CSA  open/close  solution  values 
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from  the  sequential  heuristic  as  a  partial,  starting,  feasible  solution  yields  an  optimal 
objective  value  of  83.43  which  is  2%  from  the  optimal  integer  solution  and  the  lower 
bound. 

Another  comparison  for  a  5  time  period  problem  provides  further  justification  for 
use  of  the  sequential  heuristic.  As  previously  discussed,  the  MIP  for  5  time  periods  could 
not  be  formulated  and  solved.  However,  solving  the  LP  relaxation  of  the  5  time  period 
problem  with  MINOS  gives  a  lower  bound  of  121.09.  Then,  using  the  sequential 
heuristics  ASP  and  CSA  open/close  solutions  as  a  starting  feasible  solution  as  before, 
the  5  time  period  MIP  can  be  solved  in  64.5  CPU  seconds  with  an  objective  function 
value  of  121.25.  This  objective  function  value  is  .13%  from  the  lower  bound. 

The  aforementioned  examples  imply  that  the  sequential  heuristic's  solutions  are 
"good"  enough  for  analysis  particularly  when  coupled  with  much  faster  solutions  times 
and  the  capability  to  handle  much  larger  problem  size.  The  integer  portion  of  the  se- 
quential heuristic  gives  solutions  which  are  "close"  to  the  actual  integer  solutions  pro- 
vided by  the  MIP  and  "close"  to  the  lower  bound  given  by  solving  the  LP  relaxation.  In 
general,  the  sequential  heuristic  does  not  guarantee  optimal  solutions,  only  "good"  ones 
[Ref.  22]. 

Since  the  MIP  could  not  be  formulated  and  solved  for  the  given  scenario,  all  further 
analysis  was  pursued  using  the  sequential  heuristic.  Forecasting  30  days  of  ammunition 
is  probably  unrealistic.  However,  Corps'  planning  staffs  certainly  forecast  farther  ahead 
than  three  days,  something  on  the  order  of  two  weeks.  The  main  point  is  that  a  tool  is 
presented  for  warplanning  and  actual  logistic  implementation.  The  sequential  heuristic 
is  not  confined  to  3  time  periods.  If  3  time  periods  are  appropriate  or  under  some  cir- 
cumstance applicable,  a  more  rapid  solution  time  is  achieved  using  the  heuristic  versus 
the  MIP. 

Solutions  for  the  3  time  period  MIP  indicate  a  shortage  in  lift  capacity  at  the  CSA 
and  ASP  and  a  reluctance  to  carry  inventory  from  time  period  to  time  period.  To  obtain 
a  feasible  solutions,  CSA  lift  capacity  was  increased  to  10664  STON  per  day,  ASP  lift 
capacity  increased  to  2732  STON  per  day  and  the  lower  bound  for  inventory  at  the  ASP 
reduced  to  one  day  of  supply.  In  addition,  80%  of  the  authorized  tractor/ trailers  were 
required  as  well  as  an  80%  maintenance  availability.  These  results  are  similar  to  those 
discovered  using  the  sequential  heuristic,  so  an  extended  discussion  is  delayed. 


36 


B.     THE  SEQUENTIAL  HEURISTIC 

Following  the  doctrinal  guidelines  of  Chapter  II,  appropriate  parameters  and  scalars 
were  entered  into  the  GAMS  code  of  Appendix  A  for  30  time  periods  and  executed. 

The  sequential  heuristic  requires  21.9  CPU  seconds  to  arrive  at  optimality.  The 
ammunition  facility  location  model  is  solved  in  7.4  CPU  seconds  using  ZOOM.  MINOS 
solves  the  network  flow  model  in  14.5  CPU  seconds.  A  combined  solution  time  of  21.9 
CPU  seconds  with  less  than  2  megabytes  of  computer  memory  is  credited  to  the  "such 
that"  operator.  Constraint  and  variable  reduction  yields  1861  binary  variables,  487  con- 
straints and  a  sparsity  of  1.46%  for  the  first  module;  5215  continuous  variables,  845 
constraints  and  sparsity  of  0.57%  for  the  second  module.  This  is  significant  since  early 
attempts  at  model  formulation  would  have  generated  a  combined  total  of  over  70,000 
constraints! 

By  doctrine,  using  bypass  percentages  of  80%  from  the  CSA  to  the  ATP,  lift  ca- 
pacities of  5332  STON  per  day  for  the  CSA  and  1366  STON  per  day  for  the  ASP,  an 
inventory  goal  of  3-5  days  of  supply  at  the  ASP,  75%  for  maintenance  availability  and 
75%  of  Corps  tractor/ trailers  dedicated  to  hauling  ammunition,  no  feasible  solution  ex- 
ists. Evidently,  with  increased  consumption  rates  furnished  by  USAOMMCS  (the  re- 
sults simulating  an  Airland  Battle  scenario),  increased  demand  and  inventory  goals 
cannot  be  supported  by  the  existing  or  proposed  WADS. 

Assuming  that  inventory  can  be  sacrificed  to  meet  demand,  adjustments  to  goals 
were  analyzed.  Reducing  the  inventory  goal's  lower  bound  of  3  days  of  supply  at  the 
ASP  to  a  lower  bound  of  1  day  of  supply  decreases  the  number  of  infeasibilities  sub- 
stantially but  not  completely.  The  heuristic  indicates  that  whenever  possible  the  mini- 
mum inventory  will  be  held  reducing  tactical  costs  in  the  objective  function.  Since  the 
scenario  demand  has  abrupt  "jumps",  enough  inventory  to  meet  the  minimum  require- 
ments for  the  first  and  succeeding  days  of  high  consumption  will  be  carried  forward. 
Otherwise,  minimum  inventory  is  carried  or  the  model  is  "bleeding  off'  inventory  to  get 
to  a  minimum  level. 

By  doctrine,  one  or  more  GS  companies  operate  at  the  CSA.  So  the  next  step  toward 
achieving  a  feasible  solution  was  to  increase  lift  capacity  at  the  CSA  by  a  factor  of  two 
to  10664  STON  per  day.  This  resulted  in  another  drop  in  the  number  of  infeasible  con- 
straints. Finally,  an  increase  in  ASP  lift  capacity  to  2732  STON  per  day  gave  an  optimal 
solution. 
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The  aforementioned  adjustments  to  achieve  optimality  are  a  long  way  from  the 
"rules  of  thumb"  developed  in  past  conflicts  and  published  in  current  Army  Field  Man- 
uals. The  important  question  to  ask.  is,  "Are  simulated  consumption  rates  an  accurate 
representation  of  reality?"  If  not,  what  is?  This  is  a  question  for  the  Army  analytic 
community  and  one  of  much  controversy  as  previously  mentioned. 

When  optimal,  feasible  solutions  were  obtained  from  the  sequential  heuristic,  the 
model's  flow  and  inventory  behavior  were  examined. 

After  an  inventory  level  is  set  at  ASP  j,  the/*,  flows  build  the  inventory.  When  in- 
ventory levels  become  too  high,  for  the  following  period's  stockage  level,  the  fIJt  flows 
"bleed"  inventory  off.  The/*,  bypass  flow  satisfies  the  majority  of  the  demand  and,  if 
allowed,  will  meet  all  demand  requirements  eliminating  the  need  for/*,  and/,,  flows  and 
/,,  inventory.  Since  there  is  a  positive  cost  associated  with  holding  inventory,  the  model 
attempts  to  reduce  inventory  to  the  minimum  levels  possible.  The  inventory  goal's  lower 
bound  forces  this  minimum  level  above  zero.  Prior  to  and  during  periods  of  high  de- 
mand, limited  lift  and  limited  transportation  assets  can  force  inventory  levels  above  the 
lower  bound.  This  model  behavior  can  be  contrasted  with  a  stochastic  inventory  model 
which  would  hold  inventory  for  a  specific  purpose  such  as  to  avoid  a  stockout  situation. 

The  heuristic's  tendency  is  to  carry  the  minimum  inventory  possible.  Inventory 
levels  hover  close  to  the  lower  bound  and  only  exceed  that  lower  bound  immediately 
prior  to,  during,  or  after  a  high  demand  period.  When  excess  inventory  is  on-hand  be- 
cause the  consumption  level  has  decreased,  the  model  will  reduce  inventory  levels  by 
meeting  demand  with  the  excess  ammunition  using  the/,,  flows. 

Varying  bypass  percentages  gives  some  interesting  results.  As  more  demand  is  met 
utilizing  the  CSA-ATP  arcs,  less  capability  is  present  to  build  and  "bleed  off'  inventory. 
At  90%  bypass  flow,  the  model  tends  to  build  and  hold  inventory  at  levels  slightly  higher 
than  at  levels  determined  with  a  doctrinal  80%  level.  As  less  demand  is  met  over  the 
bypass  arcs,  the  capability  to  increase  and  draw  down  inventory  increases.  75%  bypass 
flow  will  hold  inventory  at  a  level  very  close  to  the  minimum  requirement  or  in  a  rapid 
draw  down  to  achieve  a  minimum  level.  If  the  requirement  of  inventory  goals  is  loos- 
ened, then  this  characteristic  is  not  observed;  the  model  will  not  carry  inventory.  If  an 
inventory  goal  is  set  and  adhered  too,  then  high  bypass  percentages  hold  slightly  higher 
inventory  than  lower  bypass  percentages. 

Ideally,  the  WADS  would  only  be  composed  of  CSAs  and  ATPs.  If  precise  fore- 
casting,  sufficient  transportation,   and  steady,  predictable  enemy  activity  could  be 
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guaranteed,  no  inventory  would  be  necessary.  Ammunition  would  arrive  just  in  time  for 
use.  However,  this  digression  is  not  realistic.  Ammunition  will  be  stocked  at  the  ASP 
since  unforecasted  requirements  are  an  operational  fact  and  as  previously  mentioned 
stockouts  are  fatal.  Inventory  provides  a  level  of  security  or  buffer  against  unanticipated 
requirements. 

The  WADS  requires  all  demand  at  the  ATPs  to  be  satisfied  and  inventory  goals  at 
the  ASPs  to  be  met.  If  the  majority  of  demand  is  supplied  over  the  set  of  CSA-ATP 
arcs,  a  reduced  capability  is  left  to  build  inventory  to  required  levels.  This  is  seen  during 
brief  high  demand  periods.  Insuring  sufficient  inventory  in  accordance  with  Army  doc- 
trine, requires  larger  inventories  carried  from  period  to  period  or  a  decrease  in  the 
amount  of  demand  satisfied  on  the  bypass  arcs.  In  a  period  of  high  demand,  meeting  a 
large  portion  of  the  demand  over  the  CSA-ATP  arcs  and  meeting  doctrinal  inventory 
levels  may  cause  a  shortfall  in  CSA  lift  capacity.  When  little  lift  capacity  is  available  to 
build  inventories,  larger  inventories  are  carried  forward  with  each  time  period.  Carrying 
minimum  levels  of  inventory  at  the  ASP  allows  only  limited  recourse  in  an  emergency 
situation.  So  a  trade-off  exists;  the  WADS  can  either  meet  the  majority  of  demand  over 
the  bypass  arcs  or  meet  the  doctrinal  goals  for  inventory. 

There  are  no  easy  answers.  The  model  is  feasible  with  1  day  of  supply  at  the  ASP 
and  an  80%  bypass  requirement.  This  is  certainly  not  a  great  deal  of  security  should 
inventory  levels  be  at  a  minimum.  For  the  demand  data  given,  the  only  recourse  is  to 
provide  an  increased  lift  capability  at  the  CSAs  and  ASPs  to  overcome  inventory  short- 
ages and  insure  sufficient  transportation  assets  to  move  the  ammunition  where  needed. 

C.    TRANSPORTATION 

An  analysis  of  the  transportation  calculations  indicates  a  probable  shortfall  under 
periods  of  high  demand  depending  on  tractor/ trailer  availability  and  the  amount  of 
Corps  transportation  assets  dedicated  to  hauling  ammunition  forward.  Assuming  that 
Corps  tractor  authorization  does  not  change  and  using  the  15  STON  planning  figure  for 
average  haul  weight  of  the  22.5  ton  trailers,  the  parameters  that  can  vary  are  availability 
a  and  dedicated  assets  Pammo  which  are  related  multiplicatively  in  the  right-hand-side  of 
either  formulation.  Provided  that  a  Pammo  >  .543  the  model  will  remain  feasible.  As  this 
product  approaches  .543,  the  model  carries  more  inventory.  This  translates  into,  as  less 
transportation  is  available  due  to  maintenance  requirements  or  lack  of  dedicated  assets, 
more  inventory  must  be  carried  from  time  period  to  time  period  to  meet  inventory  goals 
and  demand.  The  model  uses  the  shorter  trips  between  the  CSA  and  ASP  as  a  more 
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efficient  use  of  transportation  to  build  inventory  for  later  demands  rather  the  longer 
CSA  to  ATP  trips  which  immediately  satisfy  demand.  As  the  product  moves  away  from 
.543  less  inventory  is  carried  forward.  These  calculations  are  not  exact  since  individual 
types  of  ammunition  sometimes  "cube  out"  before  they  weigh  out.  The  reverse  is  also 
true.  This  makes  the  15  STON  average  haul  weight  suspect  in  some  situations.  However, 
15  STON  of  ammunition  per  tractor/trailer  is  the  given  planning  figure.  In  addition,  any 
number  of  vehicles  could  be  used  to  haul  ammunition  in  an  emergency.  Impressed  ve- 
hicles would  be  taken  away  from  their  primary  tasks  moving  soldiers,  supplies  other  than 
ammunition,  etc.  However,  assuming  the  authorized  transport  assets  as  outlined  in 
Chapter  II  are  on-hand,  transportation  could  be  short  under  heavy  demand  without 
proper  management  by  the  MCC. 

D.    SCENARIO  SOLUTION 

After  execution  of  the  sequential  heuristic  with  the  given  data,  the  following  rec- 
ommendations are  presented  to  the  Corps'  Commander: 

1.  Location 

Open  ASPs  at  sites  4,  5,  6,  7,  11,  and  12  (Figure  5  on  page  41). 
Open  CSAs  at  sites  1,  2,  and  4  (Figure  5  on  page  41). 

2.  Flow 

Maintain  1-5  days  of  supply  at  each  ASP. 

Each  ASP  must  be  operated  by  one  DS  Company  instead  of  the  usual  one  DS 
Company  per  two  ASPs. 

Each  CSA  must  be  operated  by  two  GS  Companies. 

Bypass  flow  should  be  75%  for  more  control  over  inventory. 

Corps'  transportation  assets  are  marginal  during  high  demand  periods.  The 
product  of  availability  and  dedicated  assets  must  be  greater  than  0.543  for  the  Corps 
WADS  to  function  properly. 

Some  ATPs  will  exceed  lift  capability  and  require  augmentation  for  sustained 
operations. 

Accurate  forecasting  by  the  DAOs  is  critical  to  successful  support  of  the  for- 
ward brigades. 

The  Corps'  Commander,  in  turn,  could  issue  the  above  guidance  to  the  planning 
staff  for  study.  The  constraints  used  to  arrive  at  an  optimal  solution  may  not  actually 
be  tight.  Some  slack  in  the  form  of  additional  assets  may  exist  which  could  loosen  the 
formulation  and  make  implementation  easier.  However,  this  would  require  another 
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Figure  5.      Scenario  Solution 
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execution  of  the  heuristic  using  the  modified  constraints  to  determine  optimality  and 
would  require  subsequent  revised  guidance  by  the  Corps'  Commander  to  the  stafT. 

E.  MODELLING  THE  WADS  WITH  GAMS 

GAMS  is  an  excellent  tool  for  prototype  development  and  solving  small  to  medium 
sized  problems  [Ref.  17].  For  this  thesis,  the  Corps  level  was  considered  without  regard 
to  echelons  above  Corps.  For  a  thorough  analysis,  optimal  location  and  flow  from  the 
Port  to  the  ATPs  should  be  considered.  A  larger  model,  in  a  sequential  heuristic  form, 
would  have  more  stages  in  the  sequential  solving  procedure  [Ref.  14  :  Ch.  18  pp.  11-12]. 
Time  period  reduction  might  be  required. 

A  great  advantage  of  GAMS  is  the  ability  to  effortlessly  formulate  thousands  of 
equations,  view  at  least  a  portion  of  the  model  generated,  and  analyze  the  detailed  sol- 
ution report  available.  Coupled  with  on-line  error  messages,  debugging  is  much  easier 
than  other  programming  languages  [Ref.  14  :  Ch.18  pp.  1-2]. 

Flexibility  for  modelling  the  WADS'  myriad  of  "rules  of  thumb"  is  essential.  GAMS 
has  the  necessary  flexibility.  No  specific  model  form  is  presumed  or  required.  GAMS 
allows  the  user  modelling  freedom  necessary  to  analyze  the  many  facets  of  the  WADS. 

F.  REMARKS 

Theoretically,  both  the  MIP  and  sequential  heuristic  can  model  the  given  scenario 
and  arrive  at  an  optimal  solution.  Model  size,  rapid  solution  times,  and  "good"  solution 
quality  favor  the  heuristic.  A  detailed  pursuit  of  the  full  scale  model,  including  the  ech- 
elons above  Corps,  is  necessary  to  conduct  a  thorough  analysis  of  the  WADS.  The 
GAMS  formulations  presented  are  prototypes  for  initial  investigations  on  which  to  build 
larger  more  complete  models  that  include  all  facets  of  the  WADS. 
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V.     CONCLUSIONS  AND  RECOMMENDATIONS 

The  following  pages  close  this  thesis  with  conclusions  and  recommendations.  The 
questions  posed  in  Chapter  II  are  answered  and  additional  remarks  presented.  Rec- 
ommendations for  expansion  and  further  research  are  discussed. 

A.     CONCLUSIONS 

The  WADS,  at  least  for  Corps  level  and  below,  is  designed  to  support  the  force  in 
the  last  war  but  is  not  adequate  for  combat  under  an  Airland  Battle  scenario  with  the 
associated  increase  in  ammunition  consumption.  An  increase  in  lift  capacity  and  a  re- 
duction in  inventory  levels  is  required  to  obtain  feasible  solutions  that  are  either  "good" 
or  optimal.  Echelons  above  the  Corps  level  provide  little  relief  to  increased  consumption 
rates  by  the  forward  brigades  when  CSA  and  ASP  lift  capacities  are  at  a  maximum  is- 
suing ammunition. 

Transportation  is  feasible  at  doctrinal  levels.  However,  inventory  goals  and  demand 
can  only  be  met  by  close  management  of  transportation  assets.  Prior  to,  during,  and 
immediately  after  high  demand  periods,  transportation  constraints  are  tight.  As  long  as 
the  product  of  tractor  availability  and  dedicated  assets  remains  above  .543  then  sufficient 
transportation  assets  are  on-hand  and  the  WADS  is  feasible.  As  transport  assets  are  re- 
duced, more  inventory  is  carried  forward  with  time.  When  more  transportation  is  made 
available,  less  inventory  is  carried  forward. 

The  sequential  heuristic  provides  a  solid  approach  to  solving  the  problem  at  hand. 
With  minimum  use  of  CPU  time,  moderate  sized  problems  can  be  solved  with  "good" 
and  possibly  optimal  solutions.  The  MIP  on  the  other  hand  does  not  appear  to  have  any 
realistic  practical  value  when  using  GAMS  with  ZOOM  as  a  solver.  If  another  model 
generator  and  solver  were  utilized,  the  MIP  might  actually  be  used  since  the  sequential 
heuristic  takes  advantage  of  a  GAMS  solving  technique  which  may  not  be  available  to 
other  software. 

Chapter  II  introduced  some  natural  questions  that  have  arisen  as  a  result  of  the 
Airland  Battle  concept.  Each  question  is  discussed  in  order  of  presentation  from  Chapter 
II: 

Question  :  Will  the  increased  demand  for  ammunition  force  changes  in  the  current 
system?  Answer  :  Yes.  The  models  presented  indicate  a  shortage  in  lift  capacities  at  the 
CSAs,  ASPs,  and  ATPs.  Current  inventory  goals  are  not  feasible  and  must  be  reduced 
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to  obtain  optimal  network  flows.  Inference  for  echelons  above  Corps'  level  is  possible. 
When  CSA  and  ASP  lift  capacities  are  entirely  devoted  to  meeting  high  demand  at  the 
ATPs,  no  new  ammunition  can  be  received  at  the  CSA  or  the  ASP  from  the  Port  or 
TSA.  This  means  that  echelons  above  Corps  can  provide  little  assistance  when  lift  ca- 
pacity is  not  sufficient  at  lower  levels. 

Question  :  Are  the  heuristics  or  "rules  of  thumb"  from  past  wars  applicable  to  fu- 
ture conflicts?  Answer  :  No.  Heuristics  and  "rules  of  thumb"  from  past  wars  are  based 
on  consumptions  rates  far  lower  than  simulated  consumption  rates  under  the  concept 
of  Airland  Battle.  It  can  be  shown  using  either  model  presented,  that  historical  con- 
sumption rates  with  doctrinal  parameters  and  scalars  are  feasible.  Airland  Battle  con- 
sumptions rates  are  infeasible.  The  system  is  currently  designed  for  past  conflicts. 

Question  :  What  stockage  policies  can  be  used  to  minimize  the  size  of  ammunition 
facilities  thereby  reducing  target  signature?  Answer  :  Size  of  ammunition  facilities  is  di- 
rectly proportional  to  the  amount  of  ammunition  on-hand;  common  sense.  By  reducing 
bypass  flows,  inventory  levels  can  be  more  tightly  controlled  and  target  signature  re- 
duced. Current  policy  calls  for  3-5  days  of  supply  at  an  ASP.  As  shown,  an  inventory 
goal  of  3-5  days  of  supply  is  not  feasible,  while  1-5  days  of  supply  is.  The  model's 
tendancy  is  to  carry  rninimum  inventory  so  by  feasibility  criteria  alone,  inventory  and 
target  signature  is  reduced  (this  is  assuming  sufficient  lift  capacity  and  transportation 
assets  are  available).  If  the  efficiency  of  the  system  can  be  improved  particularly  in  the 
areas  of  ammunition  forecasting  and  system  control,  minimal  inventory  would  be  re- 
quired. By  increasing  lift  capacities  and  allowing  for  adequate  transportation  assets  to 
move  ammunition  forward,  inventory  levels  can  be  reduced  drastically  at  the  ASP.  A 
"rough  cut"  at  specific  items  to  stock  can  be  based  on  the  20  ammunition  types  discussed 
in  Chapter  II.  Straight  percentages  of  inventory  levels  given  by  the  sequential  heuristic 
will  provide  an  approximate  stockage  policy  by  ammunition  type  to  start  with. 

Question  :      What  is  the  best  placement  for  ammunition  storage  facilities? 
Answer  :    By  using  a  facility  location  model,  optimal  or  close  to  optimal  solutions  for 
proper  placement  can  be  determined  from  a  set  of  proposed  sites.  In  essence,  facility 
location  can  be  determined  with  respect  to  demand,  distance  between  facilities,  road 
conditions,  straight  line  distance  to  the  FLOT,  etc.  subject  to  doctrinal  "rules." 

This  concludes  the  questions  from  Chapterll.  Additional  comments  follow. 

It  would  appear  that  MOADS  is  the  proper  direction  to  pursue  in  concept.  Inven- 
tory levels  are  reduced  to  1-3  days  of  supply.  This  means  that  target  signature  for  each 
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individual  ASP  is  reduced.  By  increasing  the  amount  of  stockage  at  the  CSA,  the  same 
amount  of  ammunition  is  fielded  in  the  Corps'  area  but  in  a  more  secure  position.  In- 
creasing the  number  of  ASPs  from  2  to  3  creates  additional  dispersion  which  further  re- 
duces target  signature.  A  75%  bypass  level  allows  more  control  over  building  inventory 
to  meet  established  goals. 

PLS  appears  to  also  be  a  step  in  the  right  direction.  PLS  gives  an  overall  increase 
in  lift  capacity  network-wide  since  material  handling  equipment  is  not  necessarily  re- 
quired to  upload  or  offload  vehicles  transporting  ammunition.  This  concept  could  con- 
ceivably provide  the  lift  required  at  the  CSA  and  ASP  for  feasibility  with  current 
doctrinal  parameters.  Since  the  ATP  no  longer  exists  under  this  concept,  some  of  the  lift 
assets  authorized  to  the  ATPs  might  be  redistributed  within  the  TO. 

As  mentioned  in  Chapter  IV,  transportation  assets  appear  to  be  marginal.  Required 
amounts  of  transportation  are  computed  using  the  ammunition  flows  as  a  basis.  The 
numbers  are  not  encouraging  for  high  demand  periods.  Without  sufficient  transporta- 
tion, the  WADS  will  not  function.  Depending  on  other  than  transportation  units  for 
assets  to  make  up  the  difference  in  required  transportation,  e.g.  dump  trucks  from  engi- 
neer units,  is  poor  planning.  A  realistic,  detailed  transportation  support  plan  using  au- 
thorized assets  is  necessary  to  properly  design  the  WADS  for  Airland  Battle. 

Demand  rates  still  remain  a  subject  of  much  controversy.  Most  analysts  agree  that 
historical  consumption  rates  are  much  too  low  but  what  is  the  proper  consumption  rates 
for  Airland  Battle?  The  quality  of  Airland  Battle  consumption  rates  from  simulations 
or  estimation  directly  affects  the  design  of  WADS  to  support  the  force.  Poor  estimates 
will  result  in  a  system  that  is  over  designed  and  one  which  wastes  assets  or  one  which 
cannot  support  the  units  in  the  field.  A  distribution  for  demand  would  make  an  inven- 
tory theory  approach  possible  and  provide  a  much  more  powerful  examination  of  what 
is  required  to  support  the  forward  units.  In  addition,  if  a  mathematical  programming 
procedure  is  still  desired,  there  exists  a  capacitated,  stochastic  facility  location  model 
[Ref.  23  :  p.  273]. 

The  models  presented  are  not  confined  to  application  with  ammunition.  The  concept 
of  locating  facilities  and  determining  flow  through  those  facilities  opened  is  applicable 
to  other  logistical  areas  such  as  petroleum,  water,  food,  medical,  and  maintenance  op- 
erations to  mention  a  few.  Of  course,  the  doctrinal  constraints  would  change  but  the 
basic  ideas  would  remain  the  same. 
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B.     RECOMMENDATIONS 

Further  research  in  the  area  of  ammunition  facility  location  and  ammunition  flow 
is  worthwhile.  A  new  model  generator  and  solver  designed  to  take  advantage  of  problem 
structure  could  conceivably  be  undertaken  to  analyze  the  WADS  theater-wide  instead 
of  just  Corps  level  and  below.  GAMS  prototype  models  are  possible  for  the  WADS 
utilizing  the  MOADS  and  PLS  concepts  for  Corps  level  and  below  or  for  a  network  flow 
model  from  the  Port  to  the  ATPs  without  facility  location. 

Stochastic  Inventory  theory  could  also  provide  a  nice  approach  to  analyze  the 
WADS.  If  a  distribution  for  demand  were  available  this  approach  is  very  powerful  as  it 
accounts  for  randomness  in  ammunition  consumption  which  the  deterministic  models 
presented  do  not. 

The  dynamic  facility  location  model  [Ref.  23  :  pp.  269-270]  might  be  undertaken 
since  the  forward  storage  facilities  relocate  on  a  regular  basis.  The  dynamic  facility  lo- 
cation model  makes  siting  decisions  as  a  function  of  time.  This  approach  has  merit  par- 
ticularly for  attacking  units  over  long  periods  of  time. 

Another  idea  for  WADS  analysis  might  be  to  ignore  doctrinal  constraints  and  model 
only  structural  constraints.  This  concept  provides  an  unbiased  analysis  of  the  situation 
which  is  free  of  heuristics  and  "rules  of  thumb"  from  past  conflicts.  Conceivably,  a  better 
way  to  support  the  force  might  be  found  that  would  be  compatable  with  available  assets. 
Assuming  a  close  working  relationship  with  the  ammunition  community,  this  approach 
might  iteratively  arrive  at  a  system  that  would  best  support  the  units  in  the  field. 

Finally,  a  very  powerful  technique  that  could  be  used  for  either  a  doctrinal  or 
structural  approach  is  the  use  of  Elastic  programming.  In  an  Elastic  program,  all  con- 
straints can  be  violated  at  a  cost.  This  is  certainly  a  practical  approach  and  one  of  great 
flexability  when  solving  real  world  problems.  As  mentioned  in  Chapter  IV,  additional 
assets  might  be  found  to  loosen  the  formulation  at  some  cost.  Again,  an  iterative  devel- 
opment with  the  ammunition  community  would  conceivably  arrive  at  a  system  that 
could  improve  field  performance. 
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APPENDIX  A.  MIXED  INTEGER  PROGRAM 

$TITLE  MIXED  INTEGER  PROGRAM 

*THESIS  MODEL 

*CPT  MARK  J.  CAIN  DATE:  8  MARCH  1987 

*MODEL:  AN  AMMUNITION  FACILITY  LOCATION  AND  NETWORK  FLOW  MODEL  FOR 
*       A  CORPS  IN  THE  THEATER  OF  OPERATIONS. 

SETS  I  fixed  ammunition  transfer  points  /ATP1*ATP12/ 

J  possible  ammunition  storage  point  locations  /ASP1*ASP12/ 
K  possible  corps  storage  point  locations  /CSA1*CSA4/ 
T  time  periods  /T1*T3/ 

PARAMETER  ASPFLOT(J)  direct  distance  from  asp  j  to  the  flot 


/ASPl 
ASP2 

54 
59 

ASP3 

46 

ASP4 

41 

ASP5 

45 

ASP6 

49 

ASP7 

46 

ASP8 

42 

ASP9 

38 

ASP10 

35 

ASP11 

33 

ASP12 

40/  ; 

ARAME1 

[DR  CSAFLOT(K)  direct 

distance  from 

csa  k 

to  the  flot 

/CSA1 

86 

CSA2 

95 

CSA3 

97 

CSA4 

68/  ; 

ABLE  DIST(I,J)  road  distance  from  atp  i  to  asp 

j 

ASPl  ASP2 

ASP3 

ASP4  ASP5 

ASP6  ASP7 

ASP8  i 

\SP9 

ASP10  ASP11 

ASP12 

ATP1 

25   32 

26 

25 

30 

37   37 

42 

INF 

INF 

INF 

INF 

ATP  2 

29   23 

10 

11 

16 

23   23 

28 

INF 

INF 

INF 

INF 

ATP  3 

42   36 

23 

24 

29 

36   36 

41 

INF 

INF 

INF 

INF 

ATP4 

34   28 

15 

4 

9 

16   16 

21 

INF 

INF 

INF 

INF 

ATP5 

72   66 

53 

32 

27 

26   26 

31 

41 

50 

58 

64 

ATP6 

65   59 

46 

25 

20 

19   19 

24 

34 

43 

51 

57 

ATP  7 

72   66 

53 

32 

27 

26   26 

31 

41 

50 

58 

64 

ATP  8 

50   44 

31 

10 

5 

4    4 

9 

19 

28 

36 

42 

ATP  9 

INF  INF 

INF 

INF 

66 

65   59 

54 

44 

35 

27 

33 

ATP  10 

INF   INF 

INF 

INF 

58 

57   51 

46 

36 

27 

19 

25 

ATP  11 

INF  INF 

INF 

INF 

51 

50   44 

39 

29 

20 

12 

18 

ATP  12 

INF  INF 

INF 

INF 

41 

40   34 

29 

19 

10 

2 

8 
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TABLE  LENGTH(J,K)  road  distance  from  asp  j  to  csa  k 


CSA1  CSA2  CSA3  CSA4 


ASP1 

50 

82 

112 

112 

ASP2 

44 

76 

106 

106 

ASP3 

57 

89 

119 

119 

ASP4 

60 

79 

81 

81 

ASP5 

55 

74 

76 

76 

ASP6 

48 

67 

69 

69 

ASP7 

54 

73 

75 

75 

ASP8 

59 

78 

80 

72 

ASP9 

69 

88 

90 

62 

ASP10 

78 

84 

86 

53 

ASP11 

86 

76 

78 

45 

ASP12 

84 

70 

72 

39 

TABLE  H0WFAR(I,K)  road  distance  from  atp  i  to  csa  k 


CSA1 

CSA2 

CSA3 

CSA4 

ATP1 

75 

107 

106 

106 

ATP  2 

67 

99 

92 

92 

ATP  3 

80 

112 

105 

105 

ATP4 

74 

83 

85 

85 

ATP5 

74 

93 

95 

95 

ATP  6 

67 

86 

88 

88 

ATP  7 

74 

93 

95 

95 

ATP  8 

52 

71 

73 

73 

ATP9 

117 

103 

105 

72 

ATP  10 

109 

95 

97 

64 

ATP  11 

102 

88 

90 

57 

ATP  12 

92 

78 

80 
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*PENALTY  COSTS  ARE  CALCULATED  USING  A  LINEAR  COMBINATION  OF  ROAD 
♦CHARACTERISTICS  BASED  ON  THE  FOLLOWING  VALUES:  TWO  LANE- 1.00, 
*0NE  LANE-2.00,  AND  TRAILS -4.  00 

TABLE  PNCOSTl(I.J)  penalty  cost  for  road  from  atp  i  to  asp  j 


ASP1 

ASP2 

ASP3 

ASP4 

ASP5 

ASP6 

ASP7 

ASP8 

ASP9 

ASP10 

AS1 

311 

ASP12 

ATP1 

2.00 

2.00 

2.00 

1.56 

1.47 

1.  38 

1.46 

1.52 

4.00 

4.00 

4. 

00 

4.00 

ATP2 

2.00 

2.00 

2.00 

1.00 

1.00 

1.00 

1.  13 

1.29 

4.00 

4.00 

4. 

00 

4.00 

ATP  3 

1.69 

1.64 

1.43 

1.00 

1.45 

1.36 

1.44 

1.51 

4.00 

4.00 

4. 

00 

4.00 

ATP4 

1.80 

1.  77 

1.59 

1.00 

1.00 

1.00 

1.  19 

1.38 

4.00 

4.00 

4. 

00 

4.00 

ATP5 

1.92 

1.91 

1.89 

2.  16 

2.37 

2.42 

2.54 

2.45 

2.34 

2.28 

2. 

24 

2.  13 

ATP6 

1.69 

1.66 

1.56 

1.64 

1.80 

1.75 

2.00 

2.00 

2.00 

2.00 

2. 

00 

1.89 

ATP  7 

1.72 

1.70 

1.62 

1.72 

1.85 

1.88 

2.00 

2.00 

2.00 

2.00 

2. 

00 

1.91 

ATP8 

1.60 

1.55 

1.36 

2.00 

1.20 

1.00 

2.00 

2.00 

2.00 

2.00 

2. 

00 

1.86 

ATP9 

4.00 

4.00 

4.00 

4.00 

1.76 

1.  77 

1.80 

1.78 

1.  73 

1.66 

1. 

56 

1.45 

ATP  10 

4.00 

4.00 

4.00 

4.00 

1.  72 

1.  74 

1.76 

1.74 

1.67 

1.56 

1. 

37 

1.28 

ATP  11 

4.00 

4.00 

4.00 

4.00 

1.69 

1.  70 

1.  73 

1.69 

1.59 

1.40 

1. 

00 

1.00 

ATP  12 

4.00 

4.00 

4.00 

4.00 

1.85 

1.88 

1.94 

1.93 

1.89 

1.80 

1. 

00 

1.00  ; 

TABLE  PNC0ST2(J,K)  penalty  cost  for  road  from  asp  j  to  csa  k 
CSA1  CSA2  CSA3  CSA4 
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ASP1 

1. 

82 

1. 

61 

2. 

07 

2. 

07 

ASP2 

1. 

80 

1. 

58 

2. 

08 

2. 

08 

ASP3 

1. 

84 

1. 

64 

2. 

07 

2. 

07 

ASP4 

1. 

70 

1. 

35 

2. 

21 

2. 

21 

ASP5 

1. 

76 

1. 

38 

2. 

29 

2. 

29 

ASP6 

1. 

88 

1. 

42 

2. 

42 

2. 

42 

ASP7 

1. 

83 

1. 

42 

2. 

35 

2. 

35 

ASP8 

1. 

85 

1. 

46 

2. 

33 

2. 

29 

ASP9 

1. 

87 

1. 

52 

2. 

29 

2. 

34 

ASP10 

1. 

88 

2. 

43 

3. 

21 

2. 

40 

ASP11 

1. 

90 

2. 

47 

3. 

33 

2. 

47 

ASP12 

2. 

22 

2. 

60 

3. 

53 

2. 

69 

TABLE  PNC0ST3(I,K)  penalty  cost  for  road  from  atp  i  to  csa  k 


ATP1 
ATP  2 
ATP  3 
ATP4 
ATP5 
ATP6 
ATP  7 
ATP8 
ATP  9 
ATP  10 
ATP  11 
ATP  12 


CSA1 
1.  88 
1.87 
1.  73 
1.  78 

1.  76 
1.52 
1.57 
1.38 
2.01 
2.01 
2.01 

2.  12 


CSA2 
1.  70 
1.68 
1.60 
1.  34 
1.  70 
1.51 
1.55 
1.41 
2.23 
2.25 
2.27 
2.44 


CSA3 
2.06 
2.  07 
1.93 
2.  15 
2.42 
2.30 
2.27 
2.36 
2.  88 
2.94 
3.02 
3.28 


CSA4 


06 
07 
93 
15 
42 
30 


2.27 
2.36 
2.  13 
2.  14 
2.  16 
2.40 


TABLE  DMNATP(I,T)  demand  at  atp  i  in  time  period  t 


ATP1 
ATP  2 
ATP  3 
ATP4 
ATP5 
ATP6 
ATP  7 
ATP  8 
ATP  9 
ATP  10 
ATP  11 
ATP  12 


Tl 
124 
115 
115 
111 
124 
124 
109 
487 
354 
144 
258 
490 


T2 
124 
115 
115 
111 
124 
124 
109 
487 
354 
144 
258 
490 


T3 

.643 

.574 

.575 

.621 

.651 

.651 

.500 

2.001 

1.483 

.744 

1.215 

2.258 


SCALARS  NUMDIV 
NUMASP 
NUMCSA 
CORRES 
RLTN 
CONFIG 
PORASP 
PORCSA 
PERCSA 
TOFAR 
FARENF 
DAMNFAR 


number  of  divisions  assigned  to  the  corps  /3/ 
number  of  asps  assigned  to  the  corps  /6/ 
number  of  csas  assigned  to  the  corps  /3/ 
number  of  atps  directly  serviced  by  one  asp  /2/ 
number  of  atps  indirectly  serviced  by  one  csa  /2/ 
number  of  atps  directly  serviced  by  one  csa  /4/ 
percentage  of  demand  supplied  to  atp  by  asp  /.  2/ 
percentage  of  demand  supplied  to  asp  by  csa  /.  2/ 
percentage  of  demand  supplied  to  atp  by  csa  /.  8/ 
maximum  distance  from  atp  to  asp  /30/ 
maximum  distance  from  asp  to  csa  /100/ 
maximum  distance  from  atp  to  csa  /130/ 
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TUNE  adjustable  scalar  to  tune  the  objective  function  /.  01/ 

SCALE  put  asp  and  csa  in  terms  of  atp  distance  /20/ 

ADJUST  user  adjustable  scalar  to  shape  curve  /1.00/ 

MAXCSA  maximum  lift  capacity  of  the  csa  /10.664/ 

MAXASP  maximum  lift  capacity  of  the  asp  /2. 732/ 

TRIPS  possible  round  trips  per  day  from  csa  k  to  atp  i  /3/ 

ROUND  possible  round  trips  per  day  from  csa  k  to  asp  j  /4/ 

TRUCKS  number  of  tractor  trailers  authorized  per  Corps  /300/ 

AVGHAUL  average  haul  weight  per  trailer  /. 015/ 

AVAIL  tractor  availability  on  any  given  day  /.  80/ 

AMMO  percentage  of  total  tractor  trailers  hauling  ammo  /.  80/ 

LEVEL  work  level  devote  to  issue  at  csa  k  /. 333/ 

MAXINV  maximum  inventory  of  the  asp  in  days  of  supply  /5/ 

MININV  minimum  inventory  of  the  asp  in  days  of  supply  /l/  ; 

*THE  FOLLOWING  IS  A  PENALTY  COST  BASED  ON  A  USER  ADJUSTABLE  CONVEX 
^FUNCTION.  THE  CLOSER  AN  AMMUNITION  FACILITY  IS  TO  THE  FLOT,  THE 
*HIGHER  PENALTY  COST  PAID 

PARAMETER  PROX(J)  danger  curve  for  asp  j  ; 

PROX(J)=(SCAL£/ASPFLOT(J))**ADJUST  ; 

PARAMETER  CLOSE(K)  danger  curve  for  csa  k  ; 

CLOSE(K)=(SCALE/CSAFLOT(K))**ADJUST  ; 

DISPLAY  PROX,  CLOSE  ; 

*THE  FOLLOWING  0-1  PARAMETERS  ARE  USED  FOR  "SUCH  THAT"  OPERATORS  TO 
^SCREEN  OUT  DISTANCES  NOT  IN  ACCORDANCE  WITH  ARMY  DOCTRINE 

PARAMETER  MAXDIST(I,J)  filter  for  dist  ij  matrix  ; 
MAXDIST(I,J)  $  (DIST(I,J)  LE  T0FAR)=1  ; 

PARAMETER  MAXLENG(J.K)  filter  for  length  jk  matrix  ; 
MAXLENG(J.K)  $  (LENGTH(J,K)  LE  FARENF)=1  ; 

PARAMETER  MAXFAR(I,K)  filter  for  howfar  ik  matrix  ; 
MAXFAR(I,K)  $  (H0WFAR(I,K)  LE  DAMNFAR)=1  ; 

DISPLAY  MAXDIST,MAXLENG,MAXFAR  ; 

*THE  FOLLOWING  ARE  TEN  ERROR  CHECKS  FOR  FURTHER  SCREENING 
*0F  DATA  AND  MODEL  FORMULATION. 

*THIS  INSURES  THAT  4  ATP  LOCATIONS  HAVE  BEEN  INPUT  FOR  EACH  DIVISION. 

PARAMETER  CHECKl(I)  error  check  for  atp  index; 

CHECKl(I)  $  (CARD(I)/4  NE  NUMDIV)=1; 
PARAMETER  ERR0RCNT1  error  check  one; 

ERR0RCNT1  $  (SUM( I ,CHECK1( I))  NE  0)=1; 
ABORT  $(ERR0RCNT1)  "EXECUTION  TERMINATED  DUE  TO  ATP  INDEX  ERROR"; 

*THIS  INSURES  THAT  THERE  IS  AT  LEAST  ONE  PROPOSED  ASP  LOCATION  WITHIN  THE 
^FEASIBLE  DISTANCE  TO  AN  ATP. 

PARAMETER  CHECK2(I)  error  check  for  atp  asp  distance  feasibility; 
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CHECK2(I)  $  (SUM(J,MAXDIST(I,J))  EQ  0)=1; 
PARAMETER  ERR0RCNT2  error  check  two; 

ERR0RCNT2  $  ( SUM( I ,CHECK2( I) )  NE  0)=1; 
ABORT  $(ERR0RCNT2)  "EXECUTION  TERMINATED  NO  ASP  WITHIN  TOFAR  OF  ATP"; 

*THIS  INSURES  THAT  THERE  IS  AT  LEAST  ONE  PROPOSED  CSA  LOCATION  WITHIN  THE 
'^FEASIBLE  DISTANCE  TO  AN  ATP. 

PARAMETER  CHECK3(I)  error  check  for  atp  csa  distance  feasibility; 

CHECK3(I)  $  (SUM(K,MAXFAR(I,K))  EQ  0)=1; 
PARAMETER  ERRORCNT3  error  check  three; 

ERR0RCNT3  $  ( SUM( I ,CHECK3( I) )  NE  0)=1; 
ABORT  $(ERR0RCNT3)  "EXECUTION  TERMINATED  NO  CSA  WITHIN  DAMNFAR  OF  ATP"; 

*THIS  INSURED  THAT  THERE  IS  AT  LEAST  ONE  PROPOSED  CSA  LOCATION  WITHIN  THE 
^FEASIBLE  DISTANCE  TO  AN  ASP. 

PARAMETER  CHECK4(J)  error  check  for  asp  csa  distance  feasibility; 

CHECK4(J)  $  (SUM(K,MAXLENG(J,K))  EQ  0)=1; 
PARAMETER  ERR0RCNT4  error  check  three  a  ; 

ERR0RCNT4  $  (SUM( J,CHECK4( J) )  NE  0)=1; 
ABORT  $(ERR0RCNT4)  "EXECUTION  TERMINATED  NO  CSA  WITHIN  FARENF  OF  ASP"; 

*ASP  MUST  NOT  BE  PLACED  TOO  CLOSE  TO  THE  FLOT. 

PARAMETER  CHECK4A(J)  error  check  for  asp  flot  straight  line  distance  ; 

CHECK4A(J)  $  (ASPFLOT(J)  LE  20)=1  ; 
PARAMETER  ERR0RCNT4A  error  check  four  a  ; 

ERR0RCNT4A  $  (SUM(J,  CHECK4A(J))  NE  0)=1  ; 
ABORT  $(ERR0RCNT4A)  "EXECUTION  TERMINATED  ASP  TO  CLOSE  TO  FLOT"  ; 

*CSA  MUST  NOT  BE  PLACED  TOO  CLOSE  TO  THE  FLOT. 

PARAMETER  CHECK4B(K)  error  check  for  csa  flot  straight  line  distance  ; 

CHECK4B(K)  $  (CSAFLOT(K)  LE  50)=1  ; 
PARAMETER  ERR0RCNT4B  error  check  four  b  ; 

ERR0RCNT4B  $  (SUM(K,  CHECK4B(K))  NE  0)=1  ; 
ABORT  $(ERR0RCNT4B)  "EXECUTION  TERMINATED  CSA  TO  CLOSE  TO  FLOT"  ; 

*AS  DEFINED  IN  THE  PROGRAM,  PENALTY  COSTS  WILL  VARY  BETWEEN  1.00  AND 
*4. 00.  THE  FOLLOWING  WILL  INSURE  CORRECT  COMPUTATION. 

PARAMETER  CHECK5(I,J)  error  check  for  penalty  cost  calculations; 

CHECK5(I,J)  $  (PNC0ST1(I,J)  GT  4.0  OR  PNC0ST1(I,J)  LT  1.  0)=1; 
PARAMETER  ERR0RCNT5  error  check  five; 

ERR0RCNT5  $  (SUM(( I , J) ,CHECK5( I , J))  NE  0)=1; 
ABORT  $(ERR0RCNT5)  "EXECUTION  TERMINATED  PENALTY  COST  MISCALCULATION"; 

PARAMETER  CHECK5A(J,K)  error  check  for  penalty  cost  calculations; 

CHECK5A(J,K)  $  (PNC0ST2( J,K)  GT  4.0  OR  PNC0ST2(J,K)  LT  1. 0)=1; 
PARAMETER  ERRORCNT5A  error  check  five  a; 

ERR0RCNT5A  $  (SUM( ( J,K) ,CHECK5A( J,K) )  NE  0)=1; 
ABORT  $(ERR0RCNT5A)  "EXECUTION  TERMINATED  PENALTY  COST  MISCALCULATION"; 

PARAMETER  CHECK6(I,K)  error  check  for  penalty  cost  calculations; 
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CHECK6(I,K)  $  (PNC0ST3(I,K)  GT  4.0  OR  PNC0ST3(I,K)  LT  1.  0)=1; 
PARAMETER  ERR0RCNT6  error  check  six; 

ERR0RCNT6  $  (SUM( ( I ,K) ,CHECK6( I ,K) )  NE  0)=1; 
ABORT  $(ERR0RCNT6)  "EXECUTION  TERMINATED  PENALTY  COST  MISCALCULATION" 

*THIS  ERROR  CHECK  WILL  INSURE  THAT  THE  DEMAND  DATA  IS  IN  QUANTITIES 
*THAT  THE  ATP  CAN  REASONABLE  HANDLE. 

PARAMETER  CHECK7(I,T)  error  check  for  demand  at  atp  i  ; 

CHECK7(I,T)  $  (DMNATP(I,T)  GT  2. 500)=1  ; 
PARAMETER  ERR0RCNT7  error  check  seven  ; 

ERR0RCNT7  $  (SUM((I,T),  CHECK7( I ,T) )  NE  0)=1  ; 
ABORT  $(ERRORCNT7)  "EXECUTION  TERMINATED  DEMAND  EXCEED  ATP  CAPACITY"  ; 

*THE  FOLLOWING  FIVE  PARAMETERS  ARE  USER  ADJUSTABLE  TACTICAL  COSTS 
*F0R  MOVING  AND  HOLDING  AMMUNITION 

PARAMETER  SHIP1(J,K,T)  shipping  cost  from  csa  k  to  asp  j  in  period  t  ; 
SHIP1(J,K,T)  =  1.00  ; 

PARAMETER  SHIP2(I,J,T)  shipping  cost  from  asp  j  to  atp  i  in  period  t  ; 
SHIP2(I,J,T)=1.00  ; 

PARAMETER  SHIP3(I,K,T)  shipping  cost  from  csa  k  to  atp  i  in  period  t  ; 
SHIP3(I,K,T)=1.00  ; 

PARAMETER  INV1(J,T)  inventory  cost  at  asp  j  in  period  t  ; 
INV1(J,T)=1.00  ; 

VARIABLES  X(I,J)   asp  j  services  atp  i  ( l=yes  and  0=no) 
A(J,K)   csa  k  services  asp  j  ( l=yes  and  0=no) 
Z(I,K)   csa  k  services  atp  i  ( l=yes  and  0=no) 
Y(J)     asp  located  at  site  j  ( l=yes  and  0=no) 
B(K)     csa  located  at  site  k  ( l=yes  and  0=no) 
F(J,K,T)   flow  from  csa  k  to  asp  j  in  period  t 
TH(I,K,T)  flow  from  csa  k  to  atp  i  in  period  t 
EF(I,J,T)  flow  from  asp  j  to  atp  i  in  period  t 
IASP(J,T)  inventory  at  asp  j  at  the  end  of  time  period  t 
COST     objective  variable  ; 

BINARY  VARIABLES  X,A,Z,Y,B  ; 

POSITIVE  VARIABLES  F,  TH,  EF,  IASP  ; 

EQUATIONS 

ONESITE(I)  assign  atp  to  one  asp 

SERVICE(J)  asp  services  two  atps 

B0UND1(I,J,T)  upper  bound  for  ef  ij 

VARUPBD1(I,J)  variable  upper  bound  for  x  ij 

HANNA(J)  assign  asp  to  one  csa 

SUPPORT(K)  csa  services  two  asps 

B0UND2(J,K,T)  upper  bound  for  f  jk 

VARUPBD2(J,K)  variable  upper  bound  for  a  jk 

SINGLE(I)       assign  atp  to  one  csa 
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HELP(K) 

B0UND3(I,K,T) 

VARUPBD3(I,K) 

LIMIT 
CEILING 

BALASP(J,T) 

BALATP(I.T) 

CAPCSA(K.T) 

CAPASP(J,T) 

STABLE(K) 

UPPER(J,T) 

LOWER(J.T) 

TRANS (T) 

LONGRUN 

OBJFCN 


csa  services  four  atps 
upper  bound  for  th  ik 
variable  upper  bound  for  z  ik 

number  of  asps  are  limited  by  numasp 
number  of  csas  are  limited  by  numcsa 

flow  into  asp  j  must  equal  flow  out 

flow  into  atp  i  must  equal  flow  out 

lift  capacity  of  csa 

lift  capacity  of  asp 

surrogate  long  run  flow  balance  for  csa  k 

maximum  inventory  at  asp 

minimum  inventory  at  asp 

trans  assets  available  to  haul  ammo  in  time  period  t 

throughput  long  run  contribution  to  demand 

definition  of  cost  ; 


ONESITE(I). . 

SERVICE(J).. 

B0UND1(I,J,T), 

VARUPBD1(I,J) 

HANNA(J).. 

SUPPORT(K).. 

B0UND2(J,K,T). 

VARUPBD2(J,K) 

SINGLE(I).. 

HELP(K).  . 

B0UND3(I,K,T), 

VARUPBD3(I,K) 

LIMIT.  . 

CEILING. . 

BALASP(J,T).. 

BALATP(I,T). . 


SUM(J  $  (MAXDIST(I,J)),  X(I,J))=E=1  ; 

SUM(I  $  (MAXDIST(I.J)),  X( I , J) )=E=CORRES*Y( J)  ; 
.  EF(I,J,T)=L=MAXASP*X(I,J)  ; 
$  MAXDIST(I.J)..  X(I,J)=L=Y(J)  ; 

SUM(K  $  (MAXLENG(J,K)),  A( J,K) )=E=Y( J)  ; 

SUM(J  $  (MAXLENG(J,K)),  A( J,K))=E=RLTN*B(K)  ; 
.  F(J,K,T)=L=MAXASP*A(J,K)  ; 
$  MAXLENG(J,K)..  A( J,K)=L=B(K)  ; 

SUM(K  $  MAXFAR(I,K),  Z(I,K))=E=1  ; 

SUM(I  $  MAXFAR(I,K),  Z( I ,K))=E=CONFIG*B(K)  ; 
.  TH(I,K,T)=L=MAXCSA*Z(I,K)  ; 
$  MAXFAR(I,K)..  Z(I,K)=L=B(K)  ; 

SUM(J,Y(J))=E=NUMASP  ; 

SUM(K,B(K))=E=NUMCSA  ; 

-SUM(K  $  MAXLENG(J.K),  F( J,K,T))+SUM(I  $  MAXDIST( I, J) , 
EF(I,J,T))-IASP(J,T-1)+IASP(J,T)=E=0  ; 

-SUM(J  $  MAXDIST(I,J),  EF( I , J,T))-SUM(K  $  MAXFAR(I,K), 
TH(I,K,T))=E=-DMNATP(I,T)  ; 


CAPCSA(K,T)..  SUM(J  $  MAXLENG( J,K) ,  F( J,K,T) )+SUM(I  $  MAXFAR(I,K), 
TH(I,K,T))=L=MAXCSA*B(K)  ; 
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CAPASP(J,T)..  SUM(K  $  MAXLENG( J,K) ,  F( J,K,T))+SUM(I  $  MAXDIST( I, J) , 
EF(I,J,T))=L=MAXASP*Y(J)  ; 

STABLE(K)..    SUM((I,T)  $  MAXFAR(I,K),  TH(I,K,T))+ 

SUM((J,T)  $  MAXLENG(J,K),  F( J,K,T) )=L=LEVEL* 
CARD(T)*MAXCSA*B(K)  ; 

UPPER(J,T)  $  (ORD(T)  LT  CARD(T) ).  .   IASP( J,T)=L=MAXINV* 

SUM(I,DMNATP(I,T+1)*X(I,J))  ; 

LOWER(J,T)  $  (ORD(T)  LT  CARD(T)). .   IASP( J,T)=G=MININV* 

SUM(I,DMNATP(I,T+1)*X(I,J))  ; 

TRANS(T)..     (SUM((I,K)  $MAXFAR( I ,K) ,  TH( I ,K,T) )/TRIPS)+ 

(SUM((J,K)  $MAXLENG(J,K),  F( J,K,T) )/ROUND)=L= 
TRUCKS* AVGHAUL* A VAIL* AMMO  ; 

LONGRUN. .      SUM((I,K,T)  $  MAXFAR(I,K),  TH( I ,K,T) )=E=PERCSA* 
SUM((I,T),  DMNATP(I,T))  ; 

OBJFCN. .        COST=E= 
SUM((J,K,T)  $  MAXLENG(J.K),  SHIP1( J,K,T)*F( J,K,T))+ 
SUM((I,J,T)  $  MAXDIST(I,J),  SHIP2( I , J,T)*EF( I , J,T) )+ 
SUM((I,K,T)  $  MAXFAR(I,K),   SHIP3( I ,K,T)*TH( I ,K,T))+ 
SUM((J,T),  INV1(J,T)*IASP(J,T))+ 
SUM(K,  CLOSE(K)*B(K))+SUM(J,  PROX( J)*Y( J) )+ 

SUM((I,J)  $  MAXDIST(I,J),  TUNE*DIST( I , J)*PNCOSTl( I, J)*X(I, J))+ 
SUM((J,K)  $  MAXLENG(J.K),  TUNE*LENGTH( J,K)*PNCOST2( J,K)*A( J,K) )+ 
SUM((I,K)  $  MAXFAR(I,K),   TUNE*HOWFAR( I,K)*PNCOST3(I,K)*Z(I,K))  ; 

MODEL  LOCATION  /ALL/  ; 

OPTION  LIMROW=12,  LIMCOL=0,  ITERLIM=10000,  OPTCR=. 10,  RESLIM=600.  00  ; 

SOLVE  LOCATION  USING  MIP  MINIMIZING  COST  ; 

DISPLAY  X.  L,  A.  L,  Z.  L,  Y.  L,  B.  L  ; 

DISPLAY  F. L,  IASP. L,  TH. L,  EF.  L  ; 

*THE  FOLLOWING  ARE  POST-OPTIMALITY  CALCULATIONS  TO  DETERMINE  THE  NUMBER 
*OF  TRACTOR  TRAILERS  REQUIRED  TO  SUPPORT  THE  NETWORK 

PARAMETER  TRANSF( J,K,T)  tractor  trailers  needed  on  jk  arc  ; 
TRANSF(J,K,T)=F.  L(J,K,T)/AVGHAUL  ; 

PARAMETER  TRANSEF(I,J,T)  tractor  trailers  needed  on  ij  arc  ; 
TRANSEF(I,J,T)=EF. L(I,J,T)/AVGHAUL  ; 

PARAMETER  TRANSTH( I,K,T)  tractor  trailers  needed  on  ik  arc  ; 
TRANSTH(I,K,T)=TH.  L(I,K,T)/AVGHAUL  ; 

PARAMETER  STRANSF(T)  total  tractor  trailers  needed  in  time  period  t  ; 
STRANSF(T)=SUM((J,K),  TRANSF( J,K,T))  ; 
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PARAMETER  STRANSEF(T)  total  tractor  trailers  needed  in  time  period  t  ; 
STRANSEF(T)=SUM((I,J),  TRANSEF( I , J,T))  ; 

PARAMETER  STRANSTH(T)  total  tractor  trailers  needed  in  time  period  t  ; 
STRANSTH(T)=SUM((I,K),  TRANSTH( I ,K,T) )  ; 

PARAMETER  TIMETOT(T)  total  corps  requirement  for  time  period  t  ; 
TIMETOT(T)=STRANSF(T)+STRANSEF(T)+STRANSTH(T)  ; 

DISPLAY  TRANSF,  STRANSF,  TRANSEF,  STRANSEF,  TRANSTH, 
STRANSTH,  TIMETOT  ; 
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APPENDIX  B.     SEQUENTIAL  HEURISTIC 

$TITLE  A  SEQUENTIAL  HEURISTIC 

*THESIS  MODEL 


*CPT  MARK  J.  CAIN 


DATE:  8  MARCH  1987 


*MODEL:  AN  AMMUNITION  FACILITY  LOCATION  AND  NETWORK  FLOW  MODEL  FOR 
*       A  CORPS  IN  THE  THEATER  OF  OPERATIONS. 

CONTEXT 

This  sequential  heuristic  models  a  Corps  level  "slice"  of  the  U.S.  Army 
Wartime  Ammunition  Distribution  System  (WADS).  The  heuristic  is  composed 
of  a  binary  integer  program  for  the  first  module  and  a  network  flow  model 
for  the  second  module.  Maximum  use  of  the  "such  that"  operator  ($)  is 
used  to  generate  only  those  variables  and  constraints  necessary  to 
properly  model  the  system.  Solution  times  and  size  of  the  problem  that 
can  be  solved  favor  this  approach.  A  mixed  integer  program  is  the  proper 
procedure  but  solution  times  in  excess  of  600  CPU  seconds  limited  practi- 
cal value  for  analysis. 
$0FFTEXT 

SETS  I  fixed  ammunition  transfer  points  /ATP1*ATP12/ 

J  possible  ammunition  storage  point  locations  /ASP1*ASP12/ 
K  possible  corps  storage  point  locations  /CSA1*CSA4/ 
T  time  periods  /T1*T30/  ; 

PARAMETER  ASPFLOT(J)  direct  distance  from  asp  j  to  the  flot 


/ASP1 

54 

ASP2 

59 

ASP3 

46 

ASP4 

41 

ASP5 

45 

ASP6 

49 

ASP7 

46 

ASP8 

42 

ASP9 

38 

ASP10 

35 

ASP11 

33 

ASP12 

40/  ; 

TY  VI  rl 

/CSA1 

H.J\J  C 

86 

CSA2 

95 

CSA3 

97 

CSA4 

68/  ; 

PARAMETER  CSAFLOT(K)  direct  distance  from  csa  k  to  the  flot 


TABLE  DIST(I,J)  road  distance  from  atp  i  to  asp  j 

ASP1  ASP2  ASP3  ASP4  ASP5  ASP6  ASP7  ASP8  ASP9  ASP10  ASP11  ASP12 
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ATP1 

25 

32 

26 

25 

30 

37 

37 

42 

INF 

INF 

INF 

INF 

ATP  2 

29 

23 

10 

11 

16 

23 

23 

28 

INF 

INF 

INF 

INF 

ATP  3 

42 

36 

23 

24 

29 

36 

36 

41 

INF 

INF 

INF 

INF 

ATP4 

34 

28 

15 

4 

9 

16 

16 

21 

INF 

INF 

INF 

INF 

ATP  5 

72 

66 

53 

32 

27 

26 

26 

31 

41 

50 

58 

64 

ATP6 

65 

59 

46 

25 

20 

19 

19 

24 

34 

43 

51 

57 

ATP  7 

72 

66 

53 

32 

27 

26 

26 

31 

41 

50 

58 

64 

ATP8 

50 

44 

31 

10 

5 

4 

4 

9 

19 

28 

36 

42 

ATP9 

INF 

INF 

INF 

INF 

66 

65 

59 

54 

44 

35 

27 

33 

ATP  10 

INF 

INF 

INF 

INF 

58 

57 

51 

46 

36 

27 

19 

25 

ATP  11 

INF 

INF 

INF 

INF 

51 

50 

44 

39 

29 

20 

12 

18 

ATP  12 

INF 

INF 

INF 

INF 

41 

40 

34 

29 

19 

10 

2 

8  ; 

TABLE  LENGTH(J,K)  road  distance  from  asp  j  to  csa  k 


CSA1 

CSA2 

CSA3 

CSA4 

ASP1 

50 

82 

112 

112 

ASP2 

44 

76 

106 

106 

ASP3 

57 

89 

119 

119 

ASP4 

60 

79 

81 

81 

ASP5 

55 

74 

76 

76 

ASP6 

48 

67 

69 

69 

ASP7 

54 

73 

75 

75 

ASP8 

59 

78 

80 

72 

ASP9 

69 

88 

90 

62 

ASP10 

78 

84 

86 

53 

ASP11 

86 

76 

78 

45 

ASP12 

84 

70 

72 

39 

TABLE  H0WFAR(I,K)  road  distance  from  atp  i  to  csa  k 


CSA1 

CSA2 

CSA3 

CSA4 

ATP1 

75 

107 

106 

106 

ATP  2 

67 

99 

92 

92 

ATP  3 

80 

112 

105 

105 

ATP4 

74 

83 

85 

85 

ATP5 

74 

93 

95 

95 

ATP  6 

67 

86 

88 

88 

ATP  7 

74 

93 

95 

95 

ATP8 

52 

71 

73 

73 

ATP9 

117 

103 

105 

72 

ATP  10 

109 

95 

97 

64 

ATP  11 

102 

88 

90 

57 

ATP  12 

92 

78 

80 

47 

^PENALTY  COSTS  ARE  CALCULATED  USING  A  LINEAR  COMBINATION  OF  ROAD 

^CHARACTERISTICS  BASED  ON  THE  FOLLOWING  VALUES:  TWO  LANE-1.00, 
*ONE  LANE-2.00,  AND  TRAILS -4. 00 

TABLE  PNC0ST1(I,J)  penalty  cost  for  road  from  atp  i  to  asp  j 

ASP1  ASP2  ASP3  ASP4  ASP5  ASP6  ASP7  ASP8  ASP9  ASP10  ASP11  ASP12 

ATP1   2.00  2.00  2.00  1.56  1.47  1.38  1.46  1.52  4.00  4.00  4.00  4.00 

ATP2   2.00  2.00  2.00  1.00  1.00  1.00  1.13  1.29  4.00  4.00  4.00  4.00 

ATP3   1.69  1.64  1.43  1.00  1.45  1.36  1.44  1.51  4.00  4.00  4.00  4.00 

ATP4   1.80  1.77  1.59  1.00  1.00  1.00  1.19  1.38  4.00  4.00  4.00  4.00 
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ATP5  1.92  1.91  1.89  2.16  2.37  2.42  2.54  2.45  2.34  2.28  2.24  2.13 
ATP6  1.69  1.66  1.56  1.64  1.80  1.75  2.00  2.00  2.00  2.00  2.00  1.89 
ATP7   1.72  1.70  1.62  1.72  1.85  1.88  2.00  2.00  2.00   2.00   2.00   1.91 


ATP8   1.60  1.55  1.36  2.00  1.20  1.00 
ATP9   4.00  4.00  4.00  4.00  1.76  1.77 


2.  00  2.  00  2.  00 
1.80  1.  78  1.  73 


ATP10  4.00  4.00  4.00  4.00  1.72  1.74  1.76  1.74  1.67 


ATP11  4.00  4.00  4.00  4.00  1.69  1.70 
ATP12  4.00  4.00  4.00  4.00  1.85  1.88 


1.  73  1.69  1.59 
1.94  1.  93  1.89 


2. 

00 

2. 

00 

1. 

86 

1. 

66 

1. 

56 

1. 

45 

1. 

56 

1. 

37 

1. 

28 

1. 

40 

1. 

00 

1. 

00 

1. 

80 

1. 

00 

1. 

00 

TABLE  PNC0ST2(J,K)  penalty  cost  for  road  from  asp  j  to  csa  k 


CSA1 

CSA2 

CSA3 

CSA4 

ASP1 

1.82 

1.61 

2.07 

2.07 

ASP2 

1.  80 

1.58 

2.08 

2.08 

ASP3 

1.  84 

1.64 

2.  07 

2.07 

ASP4 

1.  70 

1.35 

2.21 

2.21 

ASP5 

1.  76 

1.  38 

2.29 

2.  29 

ASP6 

1.  88 

1.42 

2.42 

2.42 

ASP7 

1.83 

1.42 

2.35 

2.  35 

ASP8 

1.85 

1.46 

2.33 

2.  29 

ASP9 

1.87 

1.52 

2.29 

2.34 

ASP10 

1.  88 

2.43 

3.  21 

2.40 

ASP11 

1.90 

2.47 

3.33 

2.47 

ASP12 

2.22 

2.60 

3.53 

2.69 

TABLE  PNC0ST3(I,K)  penalty  cost  for  road  from  atp  i  to  csa  k 


1 

:sai 

CSA2  CSA3  CSA4 

ATP1 

1.88 

1.  70  2. 

06  2. 

06 

ATP  2 

1.87 

1.68  2. 

07  2. 

07 

ATP  3 

1.  73 

1.60  1. 

93  1. 

93 

ATP4 

1.  78 

1.34  2. 

15  2. 

15 

ATP5 

1.76 

1.  70  2. 

42  2. 

42 

ATP  6 

1.52 

1.51  2. 

30  2. 

30 

ATP  7 

1.57 

1.55  2. 

27  2. 

27 

ATP8 

1.38 

1.41  2. 

36  2. 

36 

ATP9 

2.01 

2.23  2. 

88  2. 

13 

ATP  10 

2.01 

2.25  2. 

94  2. 

14 

ATP  11 

2.01 

2.27  3. 

02  2. 

16 

ATP  12 

2.  12 

2.44  3. 

28  2. 

40  ; 

• 

ABLE  DMNATP(I.T)  demand  at  atp  i  in 

time 

period 

t 

Tl 

T2 

T3 

T4 

T5 

T6 

T7 

T8 

T9 

T10 

Til 

ATP1 

124 

.  124 

.643 

.564 

.  197 

197 

.  197 

.  197 

.643 

126 

126 

ATP  2 

115 

115 

.574 

.503 

.  176 

176 

.  176 

.  176 

.574 

118 

118 

ATP  3 

115 

115 

.575 

.503 

.  176 

176 

.  176 

176 

.575 

118 

118 

ATP4 

111 

111 

.621 

.551 

.  193 

193 

.  193 

.  193 

.621 

113 

113 

ATP5 

124 

124 

.651 

.570 

.200 

200 

.200 

.200 

.651 

127 

127 

ATP  6 

124 

124 

.651 

.570 

.  199 

199 

.  199 

.  199 

.651 

127 

127 

ATP  7 

109 

109 

.500 

.457 

.  160 

160 

.  160 

,  160 

.500 

111 

111 

ATP  8 

487 

487  2.001 

1.936 

.678 

678 

.678 

.678 

2.001 

498 

498 

ATP9 

354 

354  : 

L.483 

1.395 

.488 

488 

.488 

.488 

1.483 

354 

354 

ATP  10 

144 

144 

.  744 

.648 

.227 

227 

.227 

.227 

.  744 

144 

144 

ATP  11 

258 

258  : 

L.215 

1.085 

.380 

380 

.  380 

.380 

1.215 

258 

258 

ATP  12 

490 

490  2.258 

2.088 

.731 

731 

.  731 

.731 

2.258 

490 

490 
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+ 

ATP1 
ATP2 
ATP  3 
ATP4 
ATP5 
ATP  6 
ATP  7 
ATP  8 
ATP  9 
ATP  10 
ATP  11 
ATP  12 


T12 
126 
118 
118 
113 
127 
127 
111 
498 
354 
144 
285 
490 


T13 
128 
119 
119 
114 
129 
129 
112 
504 
354 
144 
258 
490 


T14 
128 
119 
119 
114 
129 
129 
112 
504 
354 
144 
258 
490 


T15 
.  128 
.  119 
.  119 
.  114 
.  129 
.  129 
.  112 
.  504 
.  354 
.  144 
.  258 
.490 


T16 

T17 

T18 

.  128 

128 

128   . 

.  119 

119 

119   . 

.  119 

119 

119   . 

.  114 

114 

114   . 

.  129 

129 

129   . 

.  129 

129 

129   . 

.  112 

112 

112   . 

.504 

504 

540   . 

.  354 

354 

354   . 

.  144 

144 

144   . 

.  258 

258 

258   . 

.490 

490 

490   . 

T19 
128 
119 
119 
114 
129 
129 
112 
504 
354 
114 
258 
490 


T20 

.501 

.445 

.445 

.489 

.508 

.508 

.382 

1.489 

1.  106 

.574 

.926 

1.698 


T21 
428 
380 
380 
430 
433 
433 
345 
482 
040 
486 
810 
616 


T22 

.  428 

.  380 

.  380 

.  430 

.433 

.433 

.  345 

1.  482 

1.  040 

.486 

.  810 

1.616 


+ 

ATP1 
ATP  2 
ATP3 
ATP4 
ATP5 
ATP6 
ATP  7 
ATP  8 
ATP  9 
ATP  10 
ATP  11 
ATP  12 


T23 

.414 

.367 

.  367 

.416 

.419 

.418 

.334 

1.432 

1.040 

.486 

.810 

1.616 


T24 

.414 

.  367 

.  367 

.416 

.419 

.418 

.334 

1.432 

1.040 

.486 

.810 

1.616 


T25 

T26 

.  124 

124  . 

.  115 

115   . 

.  115 

115   . 

.  Ill 

Ill   . 

.  124 

124   . 

.  124 

124  . 

.  109 

109   . 

.487 

487   . 

.  354 

354   . 

.  144 

144   . 

.258 

258   . 

.490 

490   . 

T27 
124 
115 
115 
111 
124 
124 
109 
487 
354 
144 
258 
490 


T28 
124 
115 
115 
111 
124 
124 
109 
487 
354 
144 
258 
490 


T29 
124 
115 
115 
111 
124 
124 
109 
487 
354 
144 
258 
490 


T30 
.  124 
.  115 
.  115 
.  Ill 
.  124 
.  124 
.  109 
.487 
.  354 
.  144 
.258 
.490  ; 


SCALARS  NUMDIV 
NUMASP 
NUMCSA 
CORRES 
RLTN 
CONFIG 
PORASP 
PORCSA 
PERCSA 
TOFAR 
FARENF 
DAMNFAR 
TUNE 
SCALE 
ADJUST 


number  of  divisions  per  corps  /3/ 

number  of  asps  assigned  to  the  corps  /6/ 

number  of  csas  assigned  to  the  corps  /3/ 

number  of  atps  directly  serviced  by  one  asp  /2/ 

number  of  atps  indirectly  serviced  by  one  csa  /4/ 

number  of  atps  directly  serviced  by  one  csa  /4/ 

percentage  of  demand  supplied  to  atp  by  asp  /. 2/ 

percentage  of  demand  supplied  to  asp  by  csa  /.  2/ 

percentage  of  demand  supplied  to  atp  by  csa  /.  8/ 

maximum  distance  from  atp  to  asp  /30/ 

maximum  distance  from  asp  to  csa  /100/ 

maximum  distance  from  atp  to  csa  /130/ 

adjustable  scalar  to  tune  the  objective  function  /.  01/ 

put  asp  and  csa  in  terms  of  atp  distance  /20/ 

user  adjustable  scalar  to  shape  curve  /1.00/  ; 


*THE  FOLLOWING  IS  A  PENALTY  COST  BASED  ON  A  USER  ADJUSTABLE  CONVEX 
♦FUNCTION.  THE  CLOSER  AN  AMMUNITION  FACILITY  IS  TO  THE  FLOT,  THE 
♦HIGHER  PENALTY  COST  PAID 

PARAMETER  PROX(J)  danger  curve  for  asp  j  ; 

PR0X(J)=( SCALE/ ASPFL0T( J)) ♦♦ADJUST  ; 

PARAMETER  CLOSE(K)  danger  curve  for  csa  k  ; 

CL0SE( K)=( SCALE/CSAFLOT( K) )♦♦ ADJUST  ; 

DISPLAY  PROX,  CLOSE  ; 
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*THE  FOLLOWING  0-1  PARAMETERS  ARE  USED  FOR  "SUCH  THAT"  OPERATORS  TO 
^SCREEN  OUT  DISTANCES  NOT  IN  ACCORDANCE  WITH  ARMY  DOCTRINE 

PARAMETER  MAXDIST(I,J)  filter  for  dist  ij  matrix  ; 
MAXDIST(I,J)  $  (DIST(I,J)  LE  TOFAR)=l  ; 

PARAMETER  MAXLENG(J,K)  filter  for  length  jk  matrix  ; 
MAXL£NG(J,K)  $  (LENGTH(J,K)  LE  FARENF)=1  ; 

PARAMETER  MAXFAR(I,K)  filter  for  howfar  ik  matrix  ; 
MAXFAR(I,K)  $  (HOWFAR(I,K)  LE  DAMNFAR)=1  ; 

DISPLAY  MAXDIST,MAXLENG,MAXFAR  ; 

*THE  FOLLOWING  ARE  TEN  ERROR  CHECKS  FOR  FURTHER  SCREENING 
*OF  DATA  AND  MODEL  FORMULATION. 

*THIS  INSURES  THAT  4  ATP  LOCATIONS  HAVE  BEEN  INPUT  FOR  EACH  DIVISION. 

PARAMETER  CHECKl(I)  error  check  for  atp  index; 

CHECKl(I)  $  (CARD(I)/4  NE  NUMDIV)=1; 
PARAMETER  ERR0RCNT1  error  check  one; 

ERRORCNT1  $  ( SUM( I ,CHECK1( I ) )  NE  0)=1; 
ABORT  $(ERR0RCNT1)  "EXECUTION  TERMINATED  DUE  TO  ATP  INDEX  ERROR"; 

*THIS  INSURES  THAT  THERE  IS  AT  LEAST  ONE  PROPOSED  ASP  LOCATION  WITHIN  THE 
^FEASIBLE  DISTANCE  TO  AN  ATP. 

PARAMETER  CHECK2(I)  error  check  for  atp  asp  distance  feasibility; 

CHECK2(I)  $  (SUM(J,MAXDIST(I,J))  EQ  0)=1; 
PARAMETER  ERR0RCNT2  error  check  two; 

ERR0RCNT2  $  (SUM( I ,CHECK2( I) )  NE  0)=1; 
ABORT  $(ERR0RCNT2)  "EXECUTION  TERMINATED  NO  ASP  WITHIN  TOFAR  OF  ATP"; 

*THIS  INSURES  THAT  THERE  IS  AT  LEAST  ONE  PROPOSED  CSA  LOCATION  WITHIN  THE 
^FEASIBLE  DISTANCE  TO  AN  ATP. 

PARAMETER  CHECK3(I)  error  check  for  atp  csa  distance  feasibility; 

CHECK3(I)  $  (SUM(K,MAXFAR(I,K))  EQ  0)=1; 
PARAMETER  ERRORCNT3  error  check  three; 

ERRORCNT3  $  (SUM(I,CHECK3( I))  NE  0)=1; 
ABORT  $(ERRORCNT3)  "EXECUTION  TERMINATED  NO  CSA  WITHIN  DAMNFAR  OF  ATP"; 

*THIS  INSURED  THAT  THERE  IS  AT  LEAST  ONE  PROPOSED  CSA  LOCATION  WITHIN  THE 
^FEASIBLE  DISTANCE  TO  AN  ASP. 

PARAMETER  CHECK4(J)  error  check  for  asp  csa  distance  feasiblilty; 

CHECK4(J)  $  (SUM(K,MAXLENG(J,K))  EQ  0)=1; 
PARAMETER  ERR0RCNT4  error  check  three  a  ; 

ERRORCNT4  $  (SUM( J,CHECK4( J))  NE  0)=1; 
ABORT  $(ERR0RCNT4)  "EXECUTION  TERMINATED  NO  CSA  WITHIN  FARENF  OF  ASP"; 

*ASP  MUST  NOT  BE  PLACED  TOO  CLOSE  TO  THE  FLOT. 

PARAMETER  CHECK4A(J)  error  check  for  asp  flot  straight  line  distance  ; 
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CHECK4A(J)  $  (ASPFLOT(J)  LE  20)=1  ; 
PARAMETER  ERR0RCNT4A  error  check  four  a  ; 

ERR0RCNT4A  $  (SUM(J,  CHECK4A(J))  NE  0)=1  ; 
ABORT  $(ERR0RCNT4A)  "EXECUTION  TERMINATED  ASP  TO  CLOSE  TO  FLOT"  ; 

*CSA  MUST  NOT  BE  PLACED  TOO  CLOSE  TO  THE  FLOT. 

PARAMETER  CHECK4B(K)  error  check  for  csa  flot  straight  line  distance  ; 

CHECK4B(K)  $  (CSAFLOT(K)  LE  50)=1  ; 
PARAMETER  ERR0RCNT4B  error  check  four  b  ; 

ERR0RCNT4B  $  (SUM(K,  CHECK4B(K))  NE  0)=1  ; 
ABORT  $(ERR0RCNT4B)  "EXECUTION  TERMINATED  CSA  TO  CLOSE  TO  FLOT"  ; 

*AS  DEFINED  IN  THE  PROGRAM,  PENALTY  COSTS  WILL  VARY  BETWEEN  1.00  AND 
*4.  00.  THE  FOLLOWING  WILL  INSURE  CORRECT  COMPUTATION. 

PARAMETER  CHECK5(I,J)  error  check  for  penalty  cost  calculations; 

:HECK5(I,J)  $  (PNCOSTl(I,J)  GT  4.0  OR  PNC0ST1(I,J)  LT  1.  0)  =  1; 
PARAMETER  ERR0RCNT5  error  check  five; 

ERR0RCNT5  $  (SUM( ( I , J) ,CHECK5( I , J))  NE  0)=1; 
ABORT  $(ERR0RCNT5)  "EXECUTION  TERMINATED  PENALTY  COST  MISCALCULATION"; 

PARAMETER  CHECK5A(J,K)  error  check  for  penalty  cost  calculations; 

CHECK5A(J,K)  $  (PNCOST2( J,K)  GT  4.0  OR  PNC0ST2(J,K)  LT  1.0)  =  1; 
PARAMETER  ERR0RCNT5A  error  check  five  a; 

ERRORCNT5A  $  ( SUM( ( J,K) ,CHECK5A( J,K) )  NE  0)=1; 
ABORT  $(ERRORCNT5A)  "EXECUTION  TERMINATED  PENALTY  COST  MISCALCULATION"; 

PARAMETER  CHECK6(I,K)  error  check  for  penalty  cost  calculations; 

CHECK6(I,K)  $  (PNC0ST3(I,K)  GT  4.0  OR  PNC0ST3(I,K)  LT  1. 0)=1; 
PARAMETER  ERRORCNT6  error  check  six; 

ERRORCNT6  $  (SUM(( I ,K) ,CHECK6( I ,K) )  NE  0)=1; 
ABORT  $(ERR0RCNT6)  "EXECUTION  TERMINATED  PENALTY  COST  MISCALCULATION"  ; 

*THIS  ERROR  CHECK  WILL  INSURE  THAT  THE  DEMAND  DATA  IS  IN  QUANTITIES 
*THAT  THE  ATP  CAN  REASONABLE  HANDLE. 

PARAMETER  CHECK7(I,T)  error  check  for  demand  at  atp  i  ; 

CHECK7(I,T)  $  (DMNATP(I.T)  GT  2. 300)=1  ; 
PARAMETER  ERRORCNT7  error  check  seven  ; 

ERR0RCNT7  $  (SUM((I,T),  CHECK7(I,T))  NE  0)=1  ; 
ABORT  $(ERR0RCNT7)  "EXECUTION  TERMINATED  DEMAND  EXCEED  ATP  CAPACITY"  ; 

VARIABLES  X(I,J,K)  csa  k  services  asp  j  which  services  atp  i 

*  ( l=yes  and  0=no) 

Z(I,K)   csa  k  services  atp  i  ( l=yes  and  0=no) 
Y(J)     asp  located  at  site  j  ( l=yes  and  0=no) 
B(K)     csa  located  at  site  k  ( l=yes  and  0=no) 
COST     objective  variable  ; 

BINARY  VARIABLES  X,Z,Y,B  ; 

EQUATIONS 

ONESITE(I)  assign  atp  to  one  asp  and  one  csa 

SERVICE(J)  asp  services  two  atps 

SUPPORT(K)  csa  services  two  asps 
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SINGLE(I)  assign  atp  to  one  csa 

HELP(K)  csa  services  four  atps 

VARUPBD1(I,J,K)  variable  upper  bound  for  x  ijk 

VARUPBD2(I,J,K)  variable  upper  bound  for  x  ijk 

VARUPBD3(I,K)  variable  upper  bound  for  z  ik 

LIMIT  number  of  asps  are  limited  by  numasp 

CEILING  number  of  csas  are  limited  by  nuracsa 

OBJFCN  definition  of  cost  ; 


ONESITE(I)..     SUM((J,K)  $  (MAXDIST( I , J)  AND  MAXLENG( J,K) ) , 
X(I,J,K))=E=1  ; 

SERVICE(J). .     SUM((I,K)  $  (MAXDIST( I , J)  AND  MAXLENG( J,K) ) , 
X(I,J,K))=E=CORRES*Y(J)  ; 

SUPPORT(K)..     SUM((I,J)  $  (MAXDIST(I.J)  AND  MAXLENG( J,K) ) , 
X(I,J,K))=E=RLTN*B(K)  ; 

SINGLE(I)..     SUM(K  $  MAXFAR(I,K),  Z(I,K))=E=1  ; 

HELP(K)..       SUM(I  $  MAXFAR(I,K),  Z( I ,K))=E=CONFIG*B(K)  ; 

VARUPBD1(I,J,K)  $  (MAXDIST(I,J)  AND  MAXLENG( J,K)).  .  X(I, J,K)=L=Y( J)  ; 

VARUPBD2(I,J,K)  $  (MAXDIST( I , J)  AND  MAXLENG( J,K)).  .  X(I,J,K)=L=B(K)  ; 

VARUPBD3(I,K)  $  MAXFAR(I,K).  .  Z(I,K)=L=B(K)  ; 

LIMIT..        SUM(J,Y(J))=E=NUMASP  ; 

CEILING..       SUM(K,B(K))=E=NUMCSA  ; 

OBJFCN. .        COST=E= 
SUM((I,J,K,T)  $  (MAXDIST(I,J)  AND  MAXLENG( J,K)) ,  ((PORASP* 
PNCOSTK I , J)*DIST( I , J)*TUNE*DMNATP( I ,T) )+(P0RCSA*PNC0ST2( J,K)* 
I£NGTH(J,K)*TUNE*DMNATP(I,T)))*X(I,J,K))  +  SUM((I,K,T)  $  MAXFAR(I.K), 
PERCSA*PNC0ST3( I ,K)*HOWFAR( I ,K)*TUNE*DMNATP( I ,T)*Z( I ,K) )  + 
SUM(J,  PRua(J)*Y(J))  +  SUM(K,  CLOSE(K)*B(K) )  ; 

MODEL  LOCATION  /ONESITE,  SERVICE,  SUPPORT,  SINGLE,  HELP,  VARUPBD1, 
VARUPBD2,  VARUPBD3,  LIMIT,  CEILING,  OBJFCN/  ; 

OPTION  LIMR0W=12,  LIMCOL=0,  ITERLIM=5000,  OPTCR=0. 01; 
^OPTION  RESLIM=0.00  ; 

SOLVE  LOCATION  USING  MIP  MINIMIZING  COST  ; 

DISPLAY  X.  L,  Z.L,  Y.  L,  B.  L  ; 

*THE  FOLLOWING  FOUR  ERROR  CHECKS  INSURE  BINARY  DATA  IS  PASSED  FROM  THE 
^FACILITY  LOCATION  SUBMODEL  TO  THE  NETWORK  FLOW  SUBMODEL.  IF  ANY  DATA 
*IS  NOT  BINARY,  THE  PROGRAM  WILL  ABORT. 

PARAMETER  CHECK8( I , J,K)  error  check  for  x  ijk  support  path  ; 
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CHECK8(I,J,K)  $  (X. L(I,J,K)  NE  1  AND  X.  L(I,J,K)  NE  0)  =  1  ; 
PARAMETER  ERR0RCNT8  error  check  eight  ; 

ERR0RCNT8  $  (SUM( ( I , J,K) ,  CHECK8( I , J,K) )  NE  0)=1  ; 
ABORT  $(ERR0RCNT8)  "EXECUTION  TERMINATED  X  IS  NOT  BINARY"  ; 

PARAMETER  CHECK9(I,K)  error  check  for  z  ik  throughput  path  ; 

CHECK9(I,K)  $  (Z.L(I.K)  NE  1  AND  Z.L(I,K)  NE  0)  =  1  ; 
PARAMETER  ERR0RCNT9  error  check  nine  ; 

ERR0RCNT9  $  (SUM((I,K),  CHECK9(I,K))  NE  0)=1  ; 
ABORT  $(ERR0RCNT9)  "EXECUTION  TERMINATED  Z  IS  NOT  BINARY"  ; 

PARAMETER  CHECKIO(J)  error  check  for  y  j  facility  opening  ; 

CHECKIO(J)  $  (Y.  L(J)  NE  1  AND  Y.  L( J)  NE  0)  =  1  ; 
PARAMETER  ERRORCNT10  error  check  ten  ; 

ERRORCNT10  $  (SUM(J,  CHECKIO(J))  NE  0)=1  ; 
ABORT  $(ERRORCNT10)  "EXECUTION  TERMINATED  Y  IS  NOT  BINARY"  ; 

PARAMETER  CHECKll(K)  error  check  for  b  k  facility  opening  ; 

CHECKll(K)  $  (B.L(K)  NE  1  AND  B.L(K)  NE  0)=1  ; 
PARAMETER  ERRORCNT11  error  check  eleven  ; 

ERRORCNT11  $  (SUM(K,  CHECKll(K))  NE  0)=1  ; 
ABORT  $(ERR0RCNT11)  "EXECUTION  TERMINATED  B  IS  NOT  BINARY"  ; 

*THE  FOLLOWING  FIVE  PARAMETERS  ARE  USER  ADJUSTABLE  TACTICAL  COSTS 
*FOR  MOVING  AND  HOLDING  AMMUNITION 

PARAMETER  SHIP1(J,K,T)  shipping  cost  from  csa  k  to  asp  j  in  period  t  ; 
SHIP1(J,K,T)  =  1.00  ; 

PARAMETER  SHIP2(I,J,T)  shipping  cost  from  asp  j  to  atp  i  in  period  t  ; 
SHIP2(I,J,T)=1.00  ; 

PARAMETER  SHIP3(I,K,T)  shipping  cost  from  csa  k  to  atp  i  in  period  t  ; 
SHIP3(I,K,T)=1.00  ; 

PARAMETER  INV1(J,T)  inventory  cost  at  asp  j  in  period  t  ; 
INV1(J,T)=1.00  ; 

*THE  FOLLOWING  FIVE  SETS  WILL  BE  USED  TO  SET  UP  "SUCH  THAT" 
^OPERATORS  BASED  ON  THE  OPTIMAL  SOLUTION  OF  THE  LOCATION 
^PROBLEM  TO  REDUCE  THE  NUMBER  OF  CONSTRAINTS  IN  THE 
^NETWORK  FLOW  MODEL  TO  THE  MINIMUM  POSSIBLE 

SETS  ASSG(J,K)  each  asp  j  supported  by  one  csa  k  ; 
ASSG(J,K)=YES  $  (SUM(I,X. L(I,J,K)))  ; 

SETS  SPPT(I,J)  each  atp  i  supported  by  one  asp  j  ; 
SPPT(I,J)=YES  $  (SUM(K,X.  L(I,J,K)))  ; 

SETS  DIRC(I,K)  each  atp  i  supported  by  one  csa  k  ; 
DIRC(I,K)=YES  $  Z.L(I,K)  ; 

SETS  ASPOPEN(J)  filter  for  network  flow  and  lift  capacities  ; 
ASPOPEN(J)=YES  $  Y. L(J)  ; 

SETS  CSAOPEN(K)  filter  for  network  flow  and  lift  capacities  ; 
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CSAOPEN(K)=YES  $  B.L(K)  ; 

*THIS  SET  ESTABLISHES  A  TIME  PERIOD  IN  WHICH  INVENTORY  IS  REQUIRED 
*TO  BE  A  BOUNDED  VARIABLE 

SETS  TWNINE(T)  time  period  one  :hrough  twenty-nine  ; 
TWNINE(T)=YES  $  (ORD(T)  LT  CARD(T))  ; 

DISPLAY  ASSG,  SPPT,  DIRC,  ASPOPEN,  CSAOPEN,  TWNINE  ; 


SCALARS  MAXCSA 
MAXASP 
TRIPS 
ROUND 
TIMES 
TRUCKS 
AVGHAUL 
AVAIL 
AMMO 
LEVEL 
MAXINV 
MININV 


maximum  lift  capacity 
maximum  lift  capacity 
possible  round  trips 
possible  round  trips 
possible  round  trips 
number  of  tractor  tra 
average  haul  weight  p 
tractor  availability 
percentage  of  total  t 
work,  level  devote  to 
maximum  inventory  of 
minimum  inventory  of 


of  the  csa  /10. 664/ 

of  the  asp  /2. 732/ 
per  day  from  csa  k  to  atp  i  /3/ 
per  day  from  csa  k  to  asp  j  /4/ 
per  day  from  asp  j  to  atp  i  /12/ 
ilers  authorized  per  Corps  /300/ 
er  trailer  /. 015/ 
on  any  given  day  /.  80/ 
ractor  trailers  hauling  ammo  /.  80/ 
issue  at  csa  over  time  /. 333/ 
the  asp  in  days  of  supply  /5/ 
the  asp  in  days  of  supply  /!/  ; 


VARIABLES  F(J,K,T)  flow  from  csa  k  to  asp  j  in  period  t 

TH(I,K,T)  flow  from  csa  k  to  atp  i  in  period  t 

EF(I,J,T)  flow  from  asp  j  to  atp  i  in  period  t 

IASP(J,T)  inventory  at  asp  j  at  the  end  of  time  period  t 

NCOST  new  objective  variable  ; 

POSITIVE  VARIABLES  F,  TH,  EF,  IASP  ; 


^ESTABLISH  INVENTORY  UPPER  AND  LOWER  BOUNDS  AS  APPROPRIATE. 
IASP. UP(J,T)  $  TWNINE(T)  =SUM( I  $  SPPT(I.J), 

DMNATP(I,T+1))*MAXINV  $  ASPOPEN(J)  ; 
IASP.L0(J,T)  $  TWNINE(T)  =SUM( I  $  SPPT(I.J), 

DMNATP(I,T+1))*MININV  $  ASPOPEN(J)  ; 


EQUATIONS 
BALASP(J.T) 
BALATP(I.T) 
CAPCSA(K,T) 
CAPASP(J.T) 
STABLE(K) 
TRANS(T) 
LONGRUN 
NOBJFCN 


flow  into  asp  j  must  equal  flow  out 

flow  into  atp  i  must  equal  flow  out 

lift  capacity  of  csa 

lift  capacity  of  asp 

surrogate  long  run  flow  balance  for  csa  k 

trans  assets  available  to  haul  ammo  in  time  period  t 

bypass  long  run  contribution  to  demand 

new  definition  of  cost  : 


BALASP(J,T)  $  ASPOPEN(J)..  -SUM(K  $  ASSG(J,K),  F( J,K,T))+SUM(I  $ 

SPPT(I,J),  EF(I,J,T))-IASP(J,T-1)+IASP(J,T)=E=0  ; 


BALATP(I,T).. 


-SUM(J  $  SPPT(I.J),  EF(I,J,T))-SUM(K  $  DIRC(I.K), 
TH(I,K,T))=E=-DMNATP(I,T)  ; 


CAPCSA(K,T)  $  CSAOPEN(K)..  SUM(J  $  ASSG(J,K),  F(J,K,T))+ 

SUM(I  $  DIRC(I.K),  TH(I,K,T))=L=MAXCSA  ; 
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CAPASP(J,T)  $  ASPOPEN(J)..  SUM(K  $  ASSG(J,K),  F(J,K,T))+ 

SUM(I  $  SPPT(I.J),  EF(I,J,T))=L=MAXASP  ; 

STABLE(K)    $  CSAOPEN(K)..  SUM((I,T)  $  DIRC(I,K),  TH(I,K,T))+ 

SUM((J,T)  $  ASSG(J,K),  F( J,K,T) )=L=LEVEL* 
CARD(T)*MAXCSA  ; 

TRANS(T)..         (SUM((I,K)  $  DIRC(I,K),  TH( I ,K,T) )/TRIPS)+ 

(SUM((J,K)  $  ASSG(J,K),  F( J ,K,T) )/ROUND)=L= 
TRUCKS* AVGHAUL* AVAIL* AMMO  ; 

LONGRUN.  .  SUM((I,K,T)  $  DIRC(I.K),  TH( I ,K,T) )=G= 

PERCSA*SUM((I,T),  DMNATP( I ,T) )  ; 

NOBJFCN.  .         NCOST=E= 

SUM((J,K,T)  $  ASSG(J,K), 

SHIP1(J,K,T)*F(J,K,T))+ 
SUM((I,J,T)  $  SPPT(I,J), 

SHIP2(I,J,T)*EF(I,J,T))+ 
SUM((I,K,T)  $  DIRC(I,K), 

SHIP3(I,K,T)*TH(I,K,T))+ 
SUM((J,T)  $  ASPOPEN(J),  INV1( J,T)*IASP( J,T) ) 


MODEL  STOCKAGE  /BALASP,  BALATP,  CAPCSA,  CAPASP, 

STABLE,  TRANS,  LONGRUN,  NOBJFCN/  ; 

OPTION  LIMR0W=4,  LIMCOL=0  ; 
*OPTION  RESLIM=0.  00  ; 

SOLVE  STOCKAGE  USING  LP  MINIMIZING  NCOST  ; 

OPTION  EJECT  ; 

DISPLAY  F.  L,  IASP.  L,  TH.  L,  EF.  L  ; 

OPTION  EJECT  ; 

*THE  FOLLOWING  ARE  POST-OPTIMALITY  CALCULATIONS  TO  DETERMINE  THE  NUMBER 
*0F  TRACTOR  TRAILERS  REQUIRED  TO  SUPPORT  THE  NETWORK 

PARAMETER  TRANSF( J,K,T)  tractor  trailers  needed  on  jk  arc  ; 
TRANSF(J,K,T)=F.  L( J, K,T)/ AVGHAUL  ; 

PARAMETER  TRANSEF( I , J,T)  tractor  trailers  needed  on  ij  arc  ; 
TRANSEF(I,J,T)=EF.  L( I, J, T)/ AVGHAUL  ; 

PARAMETER  TRANSTH(I,K,T)  tractor  trailers  needed  on  ik  arc  ; 
TRANSTH(I,K,T)=TH.  L( I ,K,T)/ AVGHAUL  ; 

PARAMETER  STRANSF(T)  total  tractor  trailers  needed  in  time  period  t  ; 
STRANSF(T)=SUM((J,K),  TRANSF( J,K,T) )  ; 

PARAMETER  STRANSEF(T)  total  tractor  trailers  needed  in  time  period  t  ; 
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STRANSEF(T)=SUM((I,J),  TRANSEF( I , J,T))  ; 

PARAMETER  STRANSTH(T)  total  tractor  trailers  needed  in  time  period  t  ; 
STRANSTH(T)=SUM((I,K),  TRANSTH( I ,K,T) )  ; 

PARAMETER  TIMETOT(T)  total  corps  requirement  for  time  period  t  ; 
TIMETOT(T)=STRANSF(T)+STRANSEF(T)+STRANSTH(T)  ; 

DISPLAY  TRANSF,  STRANSF,  TRANSEF,  STRANSEF,  TRANSTH, 
STRANSTH,  TIMETOT  ; 
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